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(_P ■ Abstract. 2-dimensional Matching Problem, which requires to find a matching of left- to right- 

Cn ' vertices in a balanced 2n- vertex bipartite graph, is a well-known polynomial problem, while various 

; , ' variants, like the 3-dimensional analogoue (3DM, with triangles on a tripartite graph), or the 

Q,' Hamiltonian Circuit Problem (HC, a restriction to "unicyclic" matchings) are among the main 

^Ji. , examples of NP-hard problems, since the first Karp reduction series of 1972. The same holds for 

^~N . the weighted variants of these problems, the Linear Assignment Problem being polynomial, and 

r^jrs ' the Numerical 3-Dimensional Matching and Travelling Salesman Problem being NP-complete. 

In this paper we show that a small modification of the 2-dimensional Matching and Assignment 
Problems in which for each i < n/2 it is required that either 7r(2j — 1) = 2« — 1 or 7r(2j) = 2i, is a 
NP-complete problem. The proof is by linear reduction from SAT (or NAE-SAT), with the size n 
of the Matching Problem being four times the number of edges in the factor graph representation 
. of the boolean problem. As a corollary, in combination with the simple linear reduction of One- 

y , J • in-Two Matching to 3-Dimensional Matching, we show that SAT can be linearly reduced to 3DM, 

2 ' while the original Karp reduction was only cubic. 
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1. Introduction 
^. 

^^ • Consider the problems encodable in tlie framework of the theory of Turing machines (say, the 

problems which could be translated into a computer program) . Say that an instance of a problem 
has size A'^ if it can be encoded into a sequence of N bits, within a given dictionary depending only 
on the problem. In a few words, Complexity Theory deals with the classification of these problems 
according to the asymptotic behaviour in N of the time of solution in the most difficult instance of 
size at most N . This kind of analysis is called worst-case analysis, and is sided by an average-case 
analysis, where the time of the solution is the average in a given measure over the ensemble of 



H \ instances of size N . 

We can restrict our attention, from this large class of problems, to problems of a special kind, 
such that, for any size N , a known finite set Sn of "feasible solutions" (or "configurations") exist, 
and, given an instance, the problem is to find the solution minimizing a certain integer-valued cost 
function (or objective function); or, if the cost function is valued in {True, False}, to find a solution 
evaluated to True or a certificate that all solutions are evaluated to False. Problems in the first and 
second class are called respectively optimization and decision problems. If, for each feasible solution, 
the time required for the evaluation of its cost is finite, any of these problems allows at least for one 
finite-time algorithm, the one which just tries sequentially all the feasible solutions. 

The class P is the class of problems such that the worst-case complexity is bounded by a given 
polynomial in A'^. Up to exotic cases, you can think to the class of problems such that, for each 
problem, a given algorithm has been explicitly provided, and has been proven to have a worst-case 
complexity bounded by the given polynomial. This class has a pleasant property: suppose that two 
problems A and B are given such that: 
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• problem B is known to be of polynomial complexity, with a certain polynomial Pb (N) ; 

• an encoding map exists, which encodes a size-iV instance of problem A into an instance of 
problem B of size at most Q{N), in a time o{Pb{Q{N))); 

• a decoding map exists, which decodes a whatever optimal solution of an instance of problem 
B, obtained from the encoding above, into an optimal solution of the original problem A, in a 
time oiPBiQiN)))- 

Then the algorithm consisting of the encoding, followed by problem- i? algorithm, and finally the 
decoding, results into a polynomial-time algorithm for problem A, with polynomial Pa{N) = 
Pb((5(7V)). This idea leads to the concept of polynomial reduction of problems. 

In the theory of Turing machines, the traditional paradigm is sided by the definition of non- 
deterministic Turing machines 0121 ■ The reader used to programming can think to these machines as 
ordinary computer programs, with the formal rule that the time of a forked process is the maximum 
of the two forked times (as if the forked processes were running on different processors), instead of 
the sum of the two times (as if the forked processes were running sequentially on a unique processor). 
The class NP of problems is the class such that the complexity, calculated within the rule above, 
is polynomial (indeed, NP stands for non- deterministic-polynomial-time). The idea of polynomial 
reduction holds also for this class of problems, as it can be applied separately to each branch of the 
process. 

Then, the crucial theorem by Cook ^3; states that any NP problem can be polynomially reduced 
to the Boolean Satisfiability Problem, i.e. to the decision problem given by n boolean variables Ui 
(i.e. literals on the alphabet {True, False}), and a statement in disjunctive normal form with m 
clauses (i.e. a string of the kind (mi V W5 V uj) A (u2 V U3 V Mg V W13) A . . . with m parentheses), a 
solution being a boolean assignment to the n variables which satisfies the statement. Despite many 
efforts, it has not been determined whether the class NP coincides with the class P (P = NP), or is 
definitely wider (PCNP), although the present intuition is maybe in the direction of stating that 
P C NP. Thus, understanding the roots of complexity commons to problems in the NP class, in order 
to compare them to the characteristics of the P class, has become an important direction of research 
of last decades. 

In particular, since every NP problem is polynomially equivalent to Boolean Satisfiability, it 
is also polynomially equivalent to any other problem which encodes Boolean Satisfiability through 
a (chain of) polynomial reductions. This induces the concept of NP-completeness for problems 
with this characteristic: via Cook theorem and the corresponding chain of reductions, a hypothetic 
polynomial-time algorithm for any NP-complete problem would automatically provide a polynomial- 
time algorithm for all problems in the NP class, thus proving that P = NP. Since the first papers by 
Karp 0|S], the family of NP-complete problems has grown, up to lead to the famous large collection 
in Garey and Johnson book jlj, and, as can be easily argued, has still grown in the 25 following 
years ^. 

1.1. Motivations and a digression 

In this paper we describe a new NP-complete problem, the One-in-Two Matching, and give a 
polynomial reduction from Boolean Satisfiability. A detailed description, section by section, of the 
contents of the paper, is postponed to the conclusions, in section|3 while here we mostly concentrate 
on the motivations. 

What is the reason, in present days, for such a work? Of course, the fact is interesting by itself, 
in the idea of making the list of NP-complctc problems still wider (a larger list of NP-complete 
problems gives a larger number of possible starting points for a reduction proof, and thus makes 
easier the task of determing whether a new problem is NP-complete). As we will see, as a side 
result we show that Onc-in-Two Matching Problem induces a chain of reductions from Boolean 
Satisfiability to the important 3-Dimensional Matching Problem (3DM in Garey- Johnson yj) which 
avoids the complicated and size-demanding original reduction by Karp g]. 

Another motivation, which was indeed the original one, is that One-in-Two Matching 
and Assignment give a hint on the structural reason why Hamiltonian Circuit (HC) and 
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Travelling Salesman (TSP) are NP-complete, although their analogue Matching and Assignment are 
polynomial, and intrinsecally simple for what concerns the energetic landscape of configurations [3 
IHlin]. Indeed, configurations of HC and TSP are permutations tt composed of a single cycle. It is 
trivial to impose that tt has no fixed points (just taking infinite weights on the diagonal), while the 
remaining restriction to have no cycles of length 2 < £ < n must be at the root of the complexity 
discrepancy among the two problems. 

In the definition of One-in-Two problems, put at infinity the non-diagonal block elements, 
W2i-i,2i = ^21.21-1 = +00, make the diagonal weight very favourite, wu — > wu — A with A — * -f oo, 
and then transpose all the row pairs (2i — 1, 2i). Then, the One-in-Two constraints is equivalent to 
forbid all the length-2 cycles among 2i — l and 2i (which are a small subset of all the cycles forbidden 
in TSP). On the other side, the matching Tr{2i) = 2z — 1, Tr{2i — 1) = 2i constitutes the obvious 
ground state of pure Assignment. 

So, although random TSP instances have good heuristics and efficient approximants, based on 
the connection with Assignment ^m]j this new ensemble of random TSP instances would be hard 
in the average case, as, instead of just condensing a few (©(Inn)) relevant long cycles, it must first 
choice how to disentangle 0{n) robust short cycles. Roughly speaking, as a cycle of length i can 
be broken in £ points, and the set of cycle lengths is a partition of n, the average complexity of this 
procedure scales with exp(n/£ ln£), and has a finite maximum for £ = 0{1). Similar arguments are 
depicted at the end of section [l~!^ , .12 , . 

Indeed, the existence of a NP-completeness proof for Onc-in-Two problems suggests that this 
narrow subset of the set of extra constraints of HC and TSP already contains the core of extra 
complexity of these problems. Also remark that, although the chain of reductions from SAT to 
Hamiltonian Circuit (3-SAT -^ Vertex-Covering -^ HC) is beautiful and elegant P, our direct 
reduction to Onc-in-Two Matching is much cheaper. 

A posteriori, maybe the main motivation for this work lies in how similar the problem is to 
situations arising in real-world experiments, to be compared to graph-, number- or set-theoretical 
problems, which require a level of mathematical abstraction. This makes One-in-Two Matching 
particularly suitable for a Quantum Adiabatic Algorithm (QAA) implementation |13lll4l[T^ll6lll7j . 

Quantum Adiabatic techniques are an alternative to more traditional Quantum Computing 
ideas. The former would solve a given NP decision problem, measuring observables on a quantum 
system which directly encodes the desired objective function (more precisely, a final-time Hamiltonian 
would reproduce the objective function, while an initial-time Hamiltonian is such that the ground 
state can be prepared, and the intermediate evolution is sufficiently slow that, by the quantum 
adiabatic theorem, with finite probability the final state encodes a solution of the problem). 

The latter, instead, would implement a quantum algorithm in a sequence of logical gates, in a 
kind of "quantum logic circuit", employing quantum bits {qubits, i.e. two-state quantum systems) 
interacting only locally, in the single gate. On one side, the existence of Quantum Error Correction 
techniques is what makes the mainstream Quantum Computing idea promising, as it would protect 
us from the unavoidable decoherence of single qubits (anyhow, cfr. J18j for Error Correction on QAA). 
The drawback is the requirement of a large number of identical quantum gates, and quantum analog 
of electric wires, to be arranged in a complex circuit, all within a support that does not add too much 
decoherence: while modern "classic" electronics allows for extremely complicated circuit patterns, 
we lack for a quantum device with a similar flexibility. 

Although the possibility that a QAA paradigm, as well as any other classical or quantum 
gedankenexperiment, could solve NP complete problems in polynomial time, is quite remote |17j . it 
has shown promising features on less ambitious tasks (cfr. mainly 16 ). So, the Quantum Adiabatic 
idea remains a promising field, provided that the mapping to the physical device is done starting 
from an "experimentally simple" NP decision problem. We can outline some points which could 
describe this vague idea: 

modularity: the easiest realizations are the ones in which it is required to design a very small 
number of elementary objects, arranged in a regular pattern. This excludes problems encoded 
on general graphs, unless we allow to deal with the (generally much larger) array which encodes 
the graph structure, i.e. for example the adjacency matrix of the graph, or the incidence matrix 
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between vertices and edges. (So, the adjacency structure of a graph, required in a graph- 
theoretical NP-complete problem, would be not cheaper than our array structure). 

low dimensionality: two-dimensional arrays of electronical devices are experimentally much more 
accessible than three-dimensional ones, and infinitely more accessible than higher-dimensional 
ones. Conversely, one-dimensional arrays, with elements having a short-range interaction, are 
probably insufficient, as the corresponding ground state is at all times a tree state [20]. Our 
array would be two-dimensional. Cfr. (1.6j for a related discussion. 

reusability: the "hard and expensive" experimental work for producing the device should allow to 
build a machine which could encode any instance of the chosen NP-complete problem, up to a 
given size. Then, there should exist a relatively "cheap and fast" machine, allowing to set up 
the specific instance, such that a given device can be used many times for different instances, 
and also for the many variants of the instance coming from the small changes of parameters 
(for example, variations in time) in the real-life problem one is interested into. In our case, the 
set of all constraints should be implemented on the device, and the set of weights should be 
encoded in tunable local "external fields" , of some physical nature, on the entries of the grid. 

small encoding: the theoretical idea of polynomial reduction could be uneffective in real-life 
applications for large polynomials (large coefficients and/or large order). One should search for 
a NP-complete problem whose reduction from the original SAT universal problem is as small 
as possible, and it is strongly adviced to be linear in the encoding size of the SAT instance, 
i.e. approximatively the number of logic operators (AND and OR) in the logic statement. We 
will see how we have a linear encoding from generic SAT, and, in section^ how the encoding 
could be still improved in real-life applications. 

no fine-tuning: the Hamiltonian should not reproduce the cost function through an exact 
cancellation of dominant terms, or a fine-tuning of contributions coming from different physical 
terms (say, like reproducing a 2-SAT clause by mean of a two-body interaction among two spins, 
and the contribution of an external magnetic field with a specific fine-tuned value which makes 
(T, T), (T, F) and {F, T) of degenerate cost), because this is unlikely to be faithfully reproduced 
by the experimental device. The unavoidable degeneracies of the cost function in a complex 
decision problem must derive uniquely from symmetry properties of the device components. 
This point will be shortly discussed in sectional 

locality: it is severely hard to reproduce non-local constraints, like the one of Hamiltonian Path or 
of Number Partitioning Problems, in a physical system, just because most of the interactions 
under experimental control have a short range. Our new One-in-Two constraint is highly local, 
as it involves occupation numbers of two {i,i) — (i + l,i -I- 1) neighbours on the grid; the 
traditional Matching constraint can be reproduced by many mechanisms, of which the most 
simple is just creating at the beginning one "row" -particle for each row (and resp. for columns), 
and implement a dynamic with conserved number of particles (cfr. section EJ . 

So, as it will be clear from the paper, contrarily to many other NP-complete problems, our Onc-in- 
Two Matching seems to meet a large number of these requirements, and it is maybe not hazardous 
to auspicate that it could be a viable candidate for a large-size experimental implementation of the 
Quantum Adiabatic Algorithm, if any will ever be pursued.:): 

1.2. Definition of 2- Dimensional Matching and Linear Assignment Problems 

Given an unoriented bipartite graph G{yi,V2;E), with Vi and V2 being the sets of the two kinds 
of vertices, both of size n, and E C Vi x V2 the set of edges, the Matching Problem j52] asks for a 
spanning subgraph G{Vi, V2; M), M C E with cardinality n, such that each vertex has degree exactly 
1, or for a certificate that such a subgraph does not exist. A set M satisfying this requirement is 

I It is worth remarking that QAA procedure on small instances has already been tried in 1191 . 



One-in-Two-Matching Problem is NP-complete 5 

called a perfect matching over G. Here and in the following, we describe a problem in the scheme 

Problem name: 



instance 
description. 

then the scheme corresponding to Matching is 

Matching Problem: 



feasible solution; 
condition to satisfy. 



GiVi,V2;E) 

with \Vi\ = \V2\ =n. 



MC E; 

\M\ ^n, V-y e ^1 U V2 degj,,f (w) = 1. 



This problem is a specific case of the more general Linear Assignment Problem, in which integer 
weights w(e) are associated to the edges, a threshold value k is given, and the search is restricted to 
perfect matchings M such that the sum of weights on the edges of M is smaller than k §. 

Linear Assignment Problem: 



G{Vi,V2;E) 
with l^il = IF2I 
w: E ^Z; 
keZ. 



MCE; 

\M\ =n, Vt; e Vi U V2 deg^f (w) 

EeGMW'(e) < k. 



More precisely, the weights could be also infinite, i.e. w : E ^ ZU {+00}, with the natural formal 
rules n + 00 = +cxd + n = +00 + 00 = +00 and +cxd > k. Then, one can assume without loss of 
generality that G is the complete balanced bipartite graph. 

In traditional notations, Matching Problem is resumed in the case w{e) = 1 for edges in E{G) 
and u'(e) = otherwise, < being replaced by >, and k = n. 

A convenient representation of these problems is via the nxn matrix W ~ {wij } of the weights. 
A feasible matching M is then described by a matrix X — {xij}, with Xij = 0, 1 and exactly one 
element equal to 1 per row and per column, such that Xij = 1 if edge e = (i, j) is in M. The cost 
function for M is restated into 



CwiX)=J2'- 



Hj •*^lj 



ivWX' 



(1) 



Another convenient representation of feasible matchings is via permutations in the symmetric group 
over n elements, tt G 5„, where 7r(i) — j if edge (i,j) is in M. In this notation the cost function 
reads 

Ci^Itt) =^w,^(j) . (2) 

i 

An example of problem instance and solution could be (on the left, items w^^u\ are written in bold) 



W = 



/3 7 2 


4 


1 


l\ 


1 6 1 


7 


8 


2 


3 3 2 


5 


6 


3 


4 2 8 


6 


2 


5 


5 5 1 


6 


3 


4 



ith k = 15; 



Mi)}i=i,....e 
C(7r) = 13. 



{3,1,2,5,6,4}; 



\4 9 8 4 3/ 



Given a whatever Assignment instance W, many algorithms allow to find in polynomial time the 
optimal assignment tt*, and its cost C* — C(7r*), for example the Hungarian Algorithm |21II22| . 



§ Equivalently, one can restrict to consider the complete balanced bipartite graph with 2n vertices, Kn,n, 
w(e) = +00 for edges in E{lCn.n) \ E{G). 



and set 
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l.S. Definition of One-in-Two Matching and One-in-Two Assignment Problems 

Now we can define the One-in-Two Matching and Assignment Problems as the variant of Matching 
(resp. Assignment) Problem in which, assumed that the dimension 2n of the matrix is even, the set 
of allowed partitions tt is restricted to include only the ones such that, for each i = 1, . . . ,n, cither 
7r(2i - 1) = 2i - 1 or 7r(2z) = 2z. 

Thus the description of One-in-Two Matching and Assignment could be resumed in the tables 



One-in-Two Matching Problem: 


G{VuV2;E) 


MCE; 


with \Vi\ = V"2 = n; 


\M\^n, Vi^GViUya deg^C^) = 1, 


partition of {Fi; V2} into 




quadruplets q = (v,v';u,u'). 


Vq {{v,u) e M) V {{v',u') e M). 



One-in-Two Assignment Problem: 



G{VuV2;E) 


MCE; 


with \Vi\ — \V2\— n; 


\M\ = n, yv cViU V2 degj^iv) = 1, 


partition of {Fi; V2} into 




quadruplets q = {v, v'; u, u'): 


Vq {{v,u)eM)\/ {{v',u')cM), 


w-.E^Z; 


J2eeMU!{e) <k. 


keZ. 





Remark that, when proven that One-in-Two Assignment is NP-complete, we will also have a proof 
that the variant with {{v,u) e M) V ((u',u') £ M) instead of {{v,u) g Af) V ((w',u') G M) 
is NP-complete. Indeed, given an instance of the OR problem, shifting the diagonal weights to 
Wii -^ Wii -\- A, and fc ^ fc + nA, we have that Cw{M) — k = n'A, with n' the number of blocks 
with two matched elements, and in the limit A — > -|-cx) we recover the analogous XOR problem. On 
the contrary, in the limit A ^ — cx) the problem becomes trivial. 

The costs of the diagonal elements do not play any role, and can be fixed to zero. Indeed, if for 
some i < n we have W2i,2i = "^21-1.21-1 — +00, no finite-cost assignment exists, while if only one of 
the two is infinite (say, W2i-i,2i-i), we are forced to fix the permutation on the other one {w2i,2i), 
and thus the elements with i or j equal to 2i never play a role: we would have had an identical cost 
function if i«2i-i,2i-i were zero, and W2i.j ~ Wj_2i — +00 for each j ^ 2i: 



W 



/ +00 


W12 


Wl3 • • • 


W21 


W22 


W23 • • • 


W31 
V : 


W32 


W33 • • • 



\ 



/ 


+ CXD 


Wl3 ... 


+ CXD 


W22 


-I-CXD 


W3I 

\ : 


-hoo 


W33 ■■■ 



(3) 



Thus, without loss of generality, we can assume that wu is finite for each i. Then, from the invariance 
of Linear Assignment, one easily convince himselfs that an equivalent instance (i.e. an instance with 
identical cost function, up to an overall constant) can be produced, with wu = for each i < 2n, 
and that the values W2i,2i-i and W2i-i,2i never appear in allowed matchings. For this reason, we 
will assume in the following that wu = 0, and denote the four elements in the n diagonal blocks of 
size 2 with special symbols * and •, instead that with a weight value. For example, an instance with 
2n — 6 could be 



W ^ 



* 


2 


4 


1 1 


* 


1 


7 


8 2 


3 3 


* 




6 3 


4 2 




* 


2 5 


5 5 


1 


6 


* 


4 9 


8 


1 


* 



with k = 10; 



The choice of representation with *s is done for mnemonic reasons: at sight, one knows that a valid 
matching should use exactly one * per block. For example, a valid matching with weight 9 could be 
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the following (on the right side, elements i such that 7r(i) = i are underlined in order to highlight 
the onc-in-two constraint satisfaction) 



/ * . 

* 


2 4 

1 7 


1 1 \ 

8 2 


3 3 

4 2 


* 


6 3 

2 5 


5 5 
\ 4 9 


1 6 

8 1 


* 



{^W}<=1...,6 = {3,2;6,4;5,1}; 
C{tt) = 9. 



Clearly, a one-in-two matching can be described by a choice of the elements kept fixed by the 
permutation (i.e. the * chosen in each block), times an allowed choice of assignment in the n- 
dimensional minor matrix resulting from the removal of the fixed rows and columns. Thus, allowed 
matchings are in bijection with pairs (it, tt), where a £ {0, 1}" and tt e iS„, with all 7r(i) ^ i. For 
example, the previous configuration could be described as (<?, {7r(i)}i=i^...^3) = ((0,0, 1), {2, 3, 1}). 
This fact suggests a naive interpretation for the potential hardness of this variant of Assignment: for 
any choice of fixed elements (the *s) , the problem of finding the optimal assignment is polynomial 
(just put +00 on the remaining *s, and use Hungarian Algorithm on the resulting Linear Assignment 
instance), nonetheless one should perform a search among these 2" possible choices, which, in 
absence of a sufficiently strong correlation or a skill mathematical structure, could make the search 
exponential in size. We will come back on this point in sectional 

2. Proof of linear reduction from Boolean Satisfiability problems 

Here we prove that also One-in-Two Matching, less general w.r.t. the analogue One-in-Two 
Assignment, allows for linear reduction from arbitrary instances of SAT, 3-SAT or NAE-3-SAT 
Problems. 

A SAT (or NAE-SAT) instance with n literals and to clauses can be encoded into a bipartite 
graph G(Vf, Vc;E), with Ve being the set of hterals {ui}i=i,...,„ and Vc the set of clauses {Ca}a=i,...,m, 
and a map s : E -^ {il} which states whether the literal enters negated or not, i.e. if Ui G Ca then 
s{i,a) = +1, while if u^ £ Ca then s{i,a) = —1. 

A satisfiability instance can be presented in the form of a factor graph, i.e. the bipartite graph 
above, where vertices in Vi are denoted by small circles, and vertices in Vc by small squares, and an 
edge is drawn in solid line if s(e) — +1, and dashed if s(e) = —1. The instance is satisfied by a given 
boolean assignment if a local constraint is satisfied on each clause vertex. For the SAT problem, 
given a clause with k incident variables, the number of solid-edge true neighbours plus dashed-edge 
false neighbours must be at least 1. For NAE-SAT it must be at least 1, and at most k — 1. 

In order to perform our reduction, it is easier to first perform a decoration on the graph: for 
each variable, introduce an auxiliary variable per incident edge (in a sense, the literal "as seen from 
the clause"), then substitute the original variable node by a "consistency check" clause (drawn as 
a small triangle), which ensures that the boolean values on the copies of the variable coincide. A 
small example of SAT factor graph, and the corresponding decorated graph, could be the following: 





In our reduction, we have two 2x2 blocks 
of the matrix W are labeled by an index (i. 



* ») for each edge (i,a) £ E{G), and thus the entries 
, a)j^2i where (1, 2) stands for the first and second block, 
and ± stands for the first and second index inside the block. Suppose to order arbitrarily the edges 
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Sa — [ion Oq), then for helping visualization, we will assume that the entries of W are ordered as 

iiii,ai)t, {ii,ai)^, ..., {ik,ak)i, {ik,ak)i, (zi,ai)^, (zi,ai)2", . . . , (ife,afe)J, {ik,ak)2) , (4) 

that is, first all the 1-blocks, then, in the same order, all the 2-blocks. 

The choice of *s in the matching corresponds to the sequence of boolean assignments for the 
literals, "as they are seen from the clause", i.e. for the literals in the decorated instance. So we 
need a "truth-setting" structure, which ensures that all these values coincide (in other words, 
implements the "consistency check" clause), and a "satisfaction-testing" structure, which checks 
that each boolean clause in the original formula is satisfied. The truth-setting structure is encoded 
in the set of Is in the entry pairs with i = j, while the satisfaction-testing structure is encoded in 
the set of Is in the entry pairs with a = b. More precisely, entries w = 1 can appear in the off-block 
matrix elements at index pairs {{i, o)„ , (j, b)g ) only in one of the two cases: 

i ^ j, a = 2 and (3=1 (truth-setting structure) 

a = b, a = 1 and (3 = 2 (satisfaction-testing structure) 
We start describing the truth-setting structures. For each variable i, call A{i) the set of adjacent 
clauses, and choose an arbitrary cyclic ordering on this set. For i ^Vi and a,b ^ A{i) we state 



y(j,a)2,(j,b)5 



1 


a = 6, 


a = T = —, 


a = 2, 


/3 = 1 


a = 6-1, 


a = T = +, 


a = 2, 


(3 = 1 





otherwise. 









(5) 



that is, the minor of W restricted to indices with fixed i looks like 



W\ 



fixed i 




/O 

/o 1 










1 0\ 

o\ 


1 




1 













1 




1 







\o 






1 



0, 

n 1/ 



(6) 



(7) 



0\ 


w' = 

*) 

Remark that in all rows with sub-index 2 and columns with sub-index 1 (i.e. all rows ands columns 
in W) we have exactly one allowed entry beyond the * element on the diagonal. In order to see 
how the truth-setting procedure works, consider what happens if we choose the top-left * in the first 
block (i.e. we choose "cr(i,ai) — 1" in the string oi(T(i_a) for the "variables seen from the clauses"). At 
the beginning the matrix is (a circle means "element chosen in the matching" , a bar means "clement 
not choosen in the matching" ) 



W 



4) 










* 






1 






* 






1 


1 






* 
* 
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while, after six logical implication we have 



W 



m 

^ 



V 



^ 



ig 



/ 



Remarking that the choice of Is in the matrix has a cyclic ordering w.r.t. clause indices, with a simple 
induction one deduces that the boolean assignments are "ci^as. — 1" for all clauses Ofe incident with 
variable i. A similar statement can be done in the case "(Ti^ai = 0" (just exchange circles with bars 
in the diagrams above). So, we have determined that only consistent boolean assignment of literals 
are allowed. 

Now we can build the structures for the satisfaction testing. We are interested in the restriction 
of W to a given fixed index a, which is of the form 



W 



h 



fixed a 



w" 



h 



(8) 



We describe the matrix W" for a clause of length k with all unncgated literals, Ca = Ui^y . 



\J Ui 



All the other cases can be trivially inferred 
neighbours of the clause a, and set 



Choose a whatever literal index (say, ii) among the 



w 



{iM)1,(].a)l 



if i = j and cr — r = — 1 ; 
if * = 3 7^ *i ^iid a = T = +1; 
\i I = ii, j ^ ii, a = +1 and r = 
li I ^ ii, j = ii, (T = — 1 and r = 

otherwise 



(9) 



that is, in an extensive representation of the matrix. 



W" = 



/ 
1 


1 



1 





1 


1 
1 








1 






1 
f 


V : : : : 





(10) 



which indeed makes the game, as one can easily check. Indeed, if all literals are negated, we are left 
with matrix minor 



W^C^.F.. 











Ik-1 



which clearly does not allow for any valid matching, while, if the first literal is true we have 

1 



W^(T,. 



■) 



4-1 



II Define the transposition matrix (T^^')ij = 1 if i = j ^ {2fc, 2k — 1}, if i = 2k and j = 2fc — 1 or if j = 2fc — 1 and 
j = 2k, and zero otherwise. I.e., matrix T^^\ acting on the left, transpose rows 2fc — 1 and 2k, while acting on the 
right transpose the corresponding columns. Then, if matrix W" encodes a clause involving {u\, . . . ,ui^, . . . ,u{), the 
matrix T^'^'W'T^^' encodes the same clause on literals (mi, . . . , m^, . . . , ui). 
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which allows at least for the diagonal matching, 7r(j) — j for all j ~ I, 
is false, but one of the others (say, the h-th) is true, we have 



, k, and if the first literal 



^(F,. 



/ 




1 






h-2 








1 





1 





V • 








h-h 



which allows at least for the matching 7r(l) — h, i:{h) — 1 and 7r(j) = j otherwise. Similarly, for a 
NAE-fc-SAT clause Ca 



W(. 



Ca = {U^, V . . 


y U^ 


j^ ) A (uii V . . . V Ui^ ) we can choose 




• 


if i = i 7^ ifc and cr = r = — 1; 
ii i = i ^ ii and a = t = +1; 




1 


ii i = ii, j ^ ii, a = +1 and r = — 1 


(^,a)l.,{3,a)l = " 


\ii ^ i\. j = i\, a = —1 and r = +1 
if i = ik, j =^ ik, cr = —1 and r = +1 






if i ^ ik, j = ik, cr — +1 and t — —1 







otherwise 



(11) 



corresponding to the matrix minor 



W"^ 



/ 
1 


1 



1 





1 


1 
1 








1 






1 
1 



1 \ 



1 



1 






\ 1 1 1 



1 
0/ 



(12) 



This completes the reduction proof for SAT and NAE-SAT (and then, in particular, for 3-SAT 
and NAE-3-SAT). Indeed, the proposed encodings of a SAT and NAE-SAT clause are a special 
case of the most general fc-literal clause, in which the number of true literals must be in the range 
{hynin, ■ ■ ■ , hmax}, (also having as a special case the l-in-3-SAT problem, which is NP-complete P), 
for which a general encoding is shown in section 0] 

It is common in reduction proofs that multiple appearence of a literal in a clause requires some 
care (truth-setting and satisfaction-testing structures could damagely interfere), and one should 
make some standard comment on the fact that this case can be excluded with small effort from any 
SAT instance. This does not happen in our case. One can understand this from the fact that the 
Is in the two structures appear in different rows and columns, and logical implications which allow 
to test the performance of the structures involve only these rows/columns. 

Also remark that the reduction is linear not only in matrix size w.r.t. the original factor-graph 
size (the dimension n of the Matching matrix is four times the number of edges in the factor graph), 
but also in the number Af of non-zero entries in the matrix, which is indeed very sparse. Each literal 
of coordination k requires a truth-setting structure with 2k entries, while each clause of length k 
requires a satisfaction-testing structure with 4A; — 3 entries (6fc — 8 for a NAE clause) , thus we have 
that, for a factor graph G{Ve, Vc] E) 

AA = 6|^(G)|-3|V;(G)| SAT problem; 

N = %\E{G)\ - 8|K(G)| NAE-SAT problem. 

Finally, we also remark that the instances of One-in-Two Matching obtained as reduction from 
Satisfiability Problems (or, as we will see in section 0] from problems with clauses having more 
generic truth tables) are of a particular kind: the Is are contained only in the top-right and bottom- 
left 2n X 2n quadrants (the matrices W' and W"). We call Bipartite One-in-Two Matching this 
specialized problem. 
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3. Proof of linear reduction from One-in- Two Matching to 3 Dimensional Matching 
We recall the definition of the S-Dimensional Matching Problem Jlj: 
3-Dimensional Matching Problem: 



GiVi,V2,V3;T) 

with T CVi XV2XV3 

and \Vi\ = IV2I = IF3I ^n. 


M C T; 

\M\ = n, 

Vu e Fi u V2 u V3 


degT(w) = 


= 1. 



The object G could be called a "tripartite hypergraph" : indeed, instead of edges, it contains hyper- 
edges with three endpoints, one for each set of vertices. Equivalently to what has been done in 
section 11.21 another representation turns out to be useful, in which the allowed elements for a 
matching are encoded in an array of zeroes and ones. Now the array is three-dimensional. Thus, 
define W = {wijk}i,j,k=i,.. .,71 such that Wijk = 1 if {i,j,k) G T, with i £ Vi, j € V2 and k G V3, 
and Wijk = otherwise. A feasible matching M is then described by an array X — {xijk}, with 
^ijk — 0, 1, such that Xijk = 1 if the triangle t = (i,j, k) is in M , and then exactly one element of 
the array is equal to 1 per i, j or k fixed, i.e. 



Vfc 



1,- 



1, 



The cost function for M is then restated into 

Cw{X) 



7 , ^ijk 
i,k 



i,3,k 



(13a) 
(13b) 
(13c) 

(14) 



and X is a valid matching if Cw{X) — n. The numerical version is defined accordingly, where the 
weights Wijk are integers, and we have a threshold value for the cost (|14|l . 

Now we describe the reduction from One-in-Two Matching to 3DM. An identical reduction 
goes from One-in-Two Assignment to Numerical 3DM. Call W — {wij} our One-in-Two Matching 
instance of dimension 2n to be encoded, with 



= < 



i = j < 2n; 

i^2h, j = 2h-l, h< n 
or i = 2/i — 1, J = 2h, h < n] 



(15) 



Wi 



(e{0,l}) otherwise. 



and W^(^) our suggested output 3DM instance, also of dimension 2n. Our formal reduction is, calling 
A{k) = {2(fc - n) - 1, 2(fc - n)} for fc == n + 1, . . . , 2n, 



w. 



(3) 
ijk 







i= j = 2k - 1, 
i=j = 2k- 

zeA(fc), 3<^A{k)- 

otherwise. 



(16) 



or, more pictorially, call e^ and Wi the vectors 



W2i-1 
W2i 



(0, 

(W2i-l,l 
(W2i,l 



i-th 

,0, 1 ,0,. 



,0); 



■ ,W2j-l,2i-2,0, 0,W;2i-l,2i+l, . 
■,W2t,2t-2 ,0, 0, W2i,2i+1 :■ 



, W2i-1^2n) 
, W2i,2n ) 



(17a) 
(17b) 
(17c) 
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then Ty(^\ written as a matrix on indices {i, k), of vectors on index j, looks like 

/ ei ... wi ... \ 

6*2 . . . W2 
1^/(3) ^ 63 . . . ^3 . . . 

6*4 ... W4 



(18) 



Indeed, remark that in the planes (i, j) for k — 1, . . . , n there are only two allowed entries, whose («, j) 
coordinates correspond to the ones of the *s in the fc-th block of the original instance. Mimicking 
the One-in-Two constraint, for each block k we are forced to choose a value Cfc £ {0, 1}, with (say) 
and 1 selecting respectively the entry with even and odd indices. Then, in all {i,j) layers with 
fc = n + 1, . . . , n there are only two non-empty i-rows, (which are empty in all the other layers with 
k = n+1, . . . , 2n). Because of the choice of the vector a, exactly one of them is now forbidden. So, 
the three-dimensional constraint of choosing one element for each index i, j and k is at this point 
reduced to a two-dimensional matching constraint, as there is a bijection between unmatched layers 
k and non-empty unmatched rows i. I.e. the 3-dimensional 2n x 2n x 2n array of equation (|18|l is 
now restricted to the n x n x n array 






W'^^'^a) 

















V 



\ 



/ 



[w. 



i^h 



^^^^2,-.(,) 



(19) 



It is easily understood that the set of forbidden indices j after the choice of vector ff and the 
disposition of the entries Wij are in agreement with the picture of section 11.31 where a vector a 
determines a n-dimensional minor of the original 2n-dimensional instance, with forbidden entries on 
the diagonal. 

As a corollary of the construction we have that the reduction is linear not only for what concerns 
the size of the array (an nxn One-in-Two instance goes into an n x n x n 3DM array) , but also on the 
number of non-zero entries in the instance ^, which is proportional to the length of the bit-encoding 
of the instance. 

Putting together this result with the one of section [21 we have in turn a linear reduction from 
SAT problems to 3DM, 

TVAjDM = 10|^(G)| - 3|Vc(G)| SAT problem; 

7V'3DM = 12|^(G)| - 8|Vc(G)| NAE-SAT problem, 

which is much more economic of the cubic one first presented in Karp seminal '72 paper 0J (see 
also 0), and maybe technically simpler. 



4. Encoding of more complex boolean patterns 

We have seen how generic boolean disjunctive- form expressions containing SAT and NAE-SAT 
clauses can be encoded into a One-in-Two Matching instance of dimension four times the number 
of edges in the factor-graph representation of the boolean formula. 

We can be more ambitious, and try to encode more complex boolean structures, in order to 
further reduce the output size of One-in-Two instances for real-life problems. For example, One- 
in- Three-SAT Problem is another NP-complete problem, where the truth tables over three literals 
only contain certain triplets of combinations (say, (TFF), [FTF] and (FFT)) 1 . A One-in-Three 
clause can clearly be encoded with five 3-SAT clauses (or three 3-SAT and one NAE-3-SAT), but 

*l More precisely, the cardinality of T for 3DM equals the one of E for One-in-Two Matching, plus the 2n 
"deterministic" entries of vectors e^, i.e. inside a factor 2 if we understand that 3DM instances having planes with 
only one valid entry can be trivially reduced in size. 
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this would require 15 (or 12) edges in the factor graph. If one could produce an encoding for One- 
in-Three clauses directly with a 6-dimensional matrix, one would save a factor 5 (or 4) in the final 
dimension of the matrix. This is the spirit in which, in section we designed explicitly compact 
matrices for NAE-SAT clauses, instead of converting them into pairs of opposite SAT clauses. 

Of course this project can not be fulfilled for arbitrary clauses of generic length: it is easily 
estimated that the number of inequivalent truth tables for length-fc clauses diverges in k much more 
severely than the number of inequivalent 2fc x 2k matrices of zeroes and ones. More precisely, 
the number of possible truth tables is 2^^ '. Then, we must consider gauge orbits (or classes, 
C) w.r.t. the symmetries of the clause, which are relabeling of literal indices (a group Sk), times 
independent negations of the literals (a group (Z2)'^). If a truth table T has a residual symmetry 
group GiT), as the orbits are a quotient, their size g{C) equals the ratio between the number of 
elements of Sk x {1^2)^ , and of G{T) 

c c ^^ ' c c 

this estimate being asymptotically very reliable ^, while the number of different matrices, also in 
the generous bound which neglects gauge orbits and trivially equivalent patterns, is at most 2^*^ . 
So, our idea of encoding in small structures more complex boolean patterns should be intended as 

(i) determine the encoding for special families of generic-length clauses, as we already did for SAT, 
NAE-SAT and truth-setting structures. 

(ii) determine the encoding for all the clause classes up to a given reasonable length. This is done 
in |Appcndix A| up to size 4. For size 5 our estimate tells us that the number of classes is higher 
than 10^, and one should be very motivated to pursue this full classification. For size 6 there 
are more than 4. • lO^'' classes, so that this classification is structurally unfeasible. 

(iii) determine the encoding for specially-important clauses (as, for example, the 5-clause used for 
summing integers). 

(iv) suggest heuristics for determining a compact encoding for new clause patterns, not discussed at 
point (jil), which could emerge in a given real-life application. 

In this section we develop the idea of point Q, for four classes of clauses of generic length k: 

2-false: clauses in which only two literal assignments for (wi, . . . , Uk) evaluate to False; 

range-T: clauses in which (mi, . . . , Uk) evaluates to True if the number of True literals is in a range 
{hm\n, ■ ■ ■ , /imax}, and variants in which each literal could enter negated; 

binary threshold: clauses in which {ui, . . . , Uk), interpreted as the binary number n{u) = 2''^^ui + 
2^~'^U2 + . . ., evaluates to True if n{u) < q, with q a threshold value; 

binary distinct (g-colouring): clauses in which (lii, . . . , Uh, w/i+i, . . . , U2h), interpreted as the set 
of two binary numbers ni{u) = 2'^~-^ui + 2'^~'^U2 + . . ■ + Uh, n2{u) = 2'^~^Uh+i + 2^~'^Uh+2 + ■ ■ ■ + 
W2/1, evaluates to True if ni{u) ^ n2{u). Allows to encode g-colouring for q — 2^, and if q is not 
a power of 2, for 2^~^ < q < 2'\ jointly with a "binary threshold" clause per colour- variable. 



Indeed, from a refined analysis 

Jitrl - _ 

2'=A:! 2**;! 



*^'^^-$r, = ^ E ( E 9(c)). (21) 

/e5fcX{Z2)^ C:feG{C) 



and, calling Rf the term in parenthesis, we have Rf < 2* for simple (ij) transpositions in Sk, Rf < 2i6 for 
permutations {ij){ki) in Sk, Rf < 22 for Z2 symmetry, and for permutations (ijk) in Sk, and so on. So, the first 
correction is of order 2 -i ~ exp(— 0.173 ■ 2 ). 
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4-1. 2-False clauses 

We expect that clauses in which almost all choices of literals are evaluated to True, or to False, are 
simpler to encode. Given a clause, call T the set of True assignments of literals, and |r| its size. The 
cases |r| = 0, 1, 2^ are trivial, and the case \T\ = 2*^ — 1 is SAT. The case \T\ = 2 corresponds either 
to the truth-setting structure (if the two True vectors are opposite on the hypercube {True, False}''), 
or otherwise to a clause which is trivial in the sense of Rule 4 of [Appendix A| and for which an 
encoding is easily deduced, combining a truth-setting structure on fc'(< k) blocks, with a bunch of 
k — k' diagonal ( g ? ) blocks. 

So, the first non-trivial class of clauses we will try to encode is the set of clauses for which 
|r| ==2*^ — 2, of which NAE-SAT is the special case in which the two False vectors are opposite. 
We call them "2-false" clauses. Without loss of generality, we can assume that the two False 
vectors are {FF ■ ■ ■ F) and {F ■ ■ ■ FT ■ ■ -T), where the number of True's is h — 1, . . . ,k, with h = k 
corresponding to NAE-SAT. Actually, h = 1 is trivial in the sense of Rule 3 of [Appendix A| For 
l</i<fc— 1, an encoding is given by the matrix below 



W 






1 


1 


1 


1 


1 


1 


u u 


■ ■ ■ 








■ ■ ■ 




1 


1 
1 


... 








... 




1 





1 
1 











1 





... 


1 





... 





1 


1 


1 











1 








1 


1 


1 


1 














1 





1 


1 



/ ) 



>h 



Indeed, for each possible choice of the literals, we have either an easy certificate of absence of any 
matching, or a matching which is valid at sight, according to the table (notation: X ■ ■ ■ X means a 
sequence of all X, while ... means arbitrary literals in the interval) 



{F---F 


F---F) 


1st row — 










{F---F 


T---T) 


1st column = 











(T ... 
i-th 
[F...T... 


... ) 
... ) 


TT{i) = i; 
^(l) = z. 


7r(z) = 1, 




<3) - 


= j o.w 


{F---F 


...TF...) 


^(l) = z. 


7r(?;) = i + 1, 


7r(i + l) = l, 


Hj) - 


= j o.w 


{F---F 


...FT...) 


7r(l) = i + l. 


7r(i + 1) = i. 


<t) = 1, 


Aj) ' 


= j o.w 



Also remark that the number of Is required in the encoding, 4fc + 2/i — 5 (< 6fc — 7), is linear in k. 



4-2. Range- T clauses 

Another family of clauses we can fully classify is what we call "range-T" clauses, i.e. fc-literal clause, 
in which the number of true literals must be in the range of values {/imin, ■ ■ • ,/imax}, for generic 
< /imin < ^max < k- Special cases are SAT (/imin = 1, hmax — k), NAE-SAT (/imin ~ 1, 
/imax — k — 1), and One-in-Three-SAT (k = 3, /imin = /imax ~ 1)- 

Call / a matrix of all J ? blocks on the diagonal, and outside. It (resp. Ip) a matrix of all 
Q 5 (resp. J g ) blocks on the diagonal, and outside; then call J a (possibly rectangular) matrix of 
all Is, and Jft (resp. Jtf) a (possibly rectangular) matrix of all g J (resp. J §) blocks. 
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Then, for given /ii.2,3, /ii + /12 + ^^3 = k, consider the matrix 

(It Jft Jft \ } ^1 
Jtf I Jft } /i2 
Jtf Jtf If / } ^3 

hi /i2 /13 

For a given assignment f of hterals, call h'l, . . . ,h'^ the number of true assignments in the three 
blocks. Also call ft," ~ hi — h[. Because of permutation invariance inside the first block, we can 
equivalently assume they are the first ones of each block, and find a, k x k reduced matrix 



//0 00\ 
J J 



Wf 



J I 





J 



J J 

\ooooo/y 



Up to a permutation both on rows and columns 



Wf = 



//0 00\ 

/ J 

J J 

J J 

J / 

\ooooo/y 



}h[ 
}h'( 



}h'i 



Then, if h'{ 



> h'2 + h'^ or /13 > /12 + h'l we have a certificate that no matching exists, because the 



rows corresponding resp. to entries h'l or h'^ have too many empty columns. Conversely, if none 
of the two facts above happens, calling h' = h'l + h'2 + h'^ and h — min(ft,", ^,3), the permutation 
7r(/i' - £) ^ h' + £ + 1, n{h' + £ + 1) ^ h' - £ ior £ ^ 0, . . . , h - 1 and tt{j) = j otherwise is a valid 
matching. 

Thus the conditions for an assignment of literals to be evaluated to True are 



h'l <h'2 + h'^ 



3 ' 



which can be restated as 



h' >hi: 



K < /12 + h'l 



h' < k- /13 



(22) 



(23) 



and, as h' is the number of True literals, we recognize the "rangc-T" constraint, with /imin = hi and 
hmux — k — hs- Remark that our choice of matrix encoding for SAT and NAE-SAT coincides with 
the general recipe described here. 



4-3. Binary threshold clauses 

A third class of clauses we encode is the "binary threshold" clause. For a set of k literals (mi, . . . ,Uk), 



denote u the binary number [ui • • • Uk] ~ 2 ui + 2 U2 + . . .. The clause evaluates to True if 
u < q, with q = [qi- ■ ■ qk\ a, threshold value with binary entries qi. The clause is trivially reduced if 
q < 2*^^^, but, as our proof of the encoding is inductive, we consider the general q < 2'^ case. 
A simple computer procedure which would determine if u < g is the following * 

for(i=l, . . ,k-l) 



> 



{ if (q[i]==l AND u[i]==0) 


Return [Yes] 


if (q[i]==0 AND u[i]==l) 


Return [No] ; 


if(q[k] >= u[k]) 


Return [Yes] 


else 


Return [No] ; 



"Rcturn[x]" means "quit the whole procedure, giving x as result". 
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We will show that our matrix encoding reproduces this code. Consider the shortcuts 

A(0)=io Ail) 



0-00 
" ~ 








1-11 



— 00 
~ 1 



B{0) 
and denote Ai = A{q. 



i?(l)=ii C(0)=io ^(1)=?,? 

Bi — B(qi), C — C{qk)- A valid encoding is given by the matrix 



(24a) 
(24b) 



(A, 
1 




Vo 





A2 

1 







Az 










Ak-1 
1 



Si \ 

B2 

B3 



Bk-i 
C J 



(25) 



where the matrix minor associated to a binary number is the one in which the odd (resp. even) rows 
and columns are deleted if the binary digit is 1 (resp. 0). 

As anticipated, we prove this inductively. First notice that C encodes the clause with fc = 1: if 
q — 1 every literal assignment is valid, while if g = only u = is accepted. This encodes the third 
"if" of the program. 

Then, in the induction assume that the minor of H25|) in which the first two rows and columns 
are removed encodes the clause on the A: — 1 digits [(72 ■ • • Qk]- Then, if qi = 1 we have 



/o 
/ 1 

1 1 
1 1 







Vo 






■• 





42 








1 


A3 •• 






1 1 



B2 
B3 



Ak-1 
1 



Bk-i 

C I 



(26) 



If u\ — 1, the first row of the reduced matrix is (1, 0, . . . , 0), so that one is forced to have 7r(l) = 1, 
and the existence of a valid matching is reduced to the constraint [u2 • • • Wfc] < \q2- ■ ■ qk\- If u\ = 0, 
the first row of the reduced matrix is (0, . . . , 0, 1), and, as, regardless to {u2, ■ ■ ■ , Uk}, we have a 
whole subdiagonal filled with Is, the permutation 7r(l) — k and 7r(i) = i — 1 for i — 2, ... ,k is a 
valid matching. This encodes the first "if" of the program. Instead, if gi = we have 



1 



1 1 
1 1 








•• 





42 








1 


^3 ■■ 








vo 



Ak-i 
1 



\ 

\ 


B2 


B3 


h-1 


C J 



(27) 



If wi = 1, the first row of the reduced matrix is (0,0, ... ,0), so that no matching is possible. If 
Ml = the first row of the reduced matrix is (1,0,..., 0), and again one is forced to have 7r(l) = 1, 
and the existence of a valid matching is reduced to the constraint [u2 ■ ■ ■ Uk] < [(/2 • ■ • Qk] ■ This 
encodes the second "if" of the program, and completes the proof. The number of Is in the matrix 
is 5fc — 4 + 2 J2i Qi (< 7^ — 4), again only linear in fc. 



4.4- Binary distinct clauses 

Finally, we will encode into a matrix with k blocks the clause which compares two binary 
numbers with h = fc/2 digits ("binary-distinct" clause). More precisely, given the literals 
(mi, . . . , Uh, Uh+i, . . . , Uk), associate the two binary numbers 



ni{u) 
n2{u) 



[Ul ■ 

Wh- 



■Uh 



■U2h\ 



2 Ul 



0/1-2 

2 U2 



-Uh, 



Uh+1 + 2 ' Uh+2 + . . • + U2h 



(28a) 
(28b) 
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and say that the clause is satisfied whenever ni(u) 7^ n2{u). This clause is an "OR" of binary 
"XOR"s, i.e. u evaluates to 

i=l,...,h 

and is encoded by the matrix (it is intended that matrix entries are where not specified) 



/ 1 1 
0011 



1 
1 



1 1 
1 1 

1 
1 



\ 1 1 



1 

1 



\ 



1 

1 



1 

1 



1011 
0111 



1 
1 



1 1 
1 1 





(29) 



00/ 
00/ 



Indeed, the only part of the matrix sensible to literal assignment is the diagonal of the top-right 
quadrant, which has Is in correspondence of the unequal-digit indices, and thus is totally empty iff 
the two binary integers are equal. In this latter case we have an easy certificate, for example the 
first h rows have complexively only h~ 1 columns not totally empty. Conversely, if the i-th element 
of the diagonal is 1 (that is, if Ui ^ Uh+i), a valid permutation is 



7r(z) = h + i 
Tr{2h) = 1 
7r(j) = j + 1 
7r(j) = j 



1 < J < «, 
i < j < h + i. 



h + i < j <2h] 



(30) 



This proves the encoding. Again remark the linear number of Is: 7k — 8. 



5. The roots of NP-hardness: a negative remark 

An astonishing fact of Complexity Theory is the frequent similarity in the pictorial description of 
pairs of problems, which are indeed one proven to be polynomial, the other to be NP-complete. 
Examples are Chinese Postman vs. Travelling Salesman, or Arc Covering vs. Vertex Covering, 
or Min-Cut-Max-Flow vs. Max-Cut-Max-Flow. In other cases, a natural threshold value seems 
to exists, for example, 2-SAT, 2-Coloring, Minimum Spanning Tree and 2-Dimcnsional Matching 
are polynomial, while K-SAT, X-Coloring, Minimum Spanning (iir-)hypertree and iiT-Diniensional 
Matching are NP-hard for any A' > 3. A generic lesson is that we lack for a reasonable intuition of 
which problems are hard, and which not. 

On the other side, a hypothetic problem of cost minimization, in which an instance consists of a 
random extraction of 2" i.i.d. "costs" variables for the 2" feasible solutions in {0, 1}" (called REM, 
Random Energy Model |23[l24ll^ ). would certainly be exponentially hard: no strategy is possible, 
behind the trivial one of just reading all the costs, and finding the smallest one, and this would take 
2" operations. Similarly, a decision problem in which each of the 2" configurations is a solution 
independently from the others, with probability p — 2^"", < a < 1, would be exponentially hard: 
no strategy is possible, behind the trivial one of just searching sequentially for a solution, the average 
number of operations required being 2(i~")". 

Of course, these problems are out of the NP class, as the mere encoding of the instance is 
exponential in size (if you want, this is the content of the celebrated Shannon entropy theorem). 
This "physical" intuition corresponds to a well-known statement of Complexity Theory, namely that 
P^NP "in an oracle setting" j'2"6j . 

The widespread idea behind (average-case) NP-hardness is that, although an ensemble of 
random instances of a given NP-complete decision problem only visits an infinitesimal fraction of 
the 2*^^"-' possible SAT-UNSAT truth tables (a subset of the ones which can be optimally encoded 
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with a polynomial entropy, as the polynonial encoding of the instance is a bound for this value), it 
is possible that they are both "sufficiently well distributed" and "sufficiently typical" , such that the 
resulting ensemble of tables is "well approximated" by the one of the REM (or "oracle" ) ensemble. 

Some interesting efforts have been recentely done in the direction of proving that certain NP- 
complete problems have statistical properties strongly similar (although not identical) to the ones of 
REM j27", 28, 29 , and this fact has been interpreted as a hint toward NP-completeness. Some other 
observations have been raised on the fact that the same kind of properties arises also in well-known 
polynomial problems, such as the determination of the rank of a matrix |8()j . 

An observation of this last nature can be done also in this case. That is, we want to state 
that not only One-in-Two Matching is formulated in an apparently similar way to other Matching 
problems which are indeed polynomial, but also the statistical properties of the spectra of costs for 
random instances of these problems share many common features. 

We know that pure Matching is polynomial. We have just proven that One-in-Two-Matching 
is NP-complete, and we gave at the end of section [T3I an argument for this fact: changing also a 
single "spin" determination ai translates into the substitution of a whole row and a whole column in 
the induced matching subproblem. What "almost always" jj happens for a random instance of, say, 
i.i.d. real positive entries Wij with measure p{w) = e"'", is that at fixed a the optimal permutation 
TT*{(j) is such that Ci4/(ct, 7r*(tT)) ~ ^(2) [^ IHll OH] and all the summands are of order l/N, the 
largest one being of order ln{N)/N. Then, taking a a' which differs from ct on a single entry, 
resamples both a whole row and a whole column. In correspondence with the old permutation we 
will now have a huge energy, (Cvi^((j',7r*((T))) ~ ({2) + 2, well inside the high-density part of the 
spectrum. Then, correlations among optimal permutations t^*{(j) and tt*{(t'), and among energies 
Cvi/(ct, 7r*((7)) — C(2) and Cw{<^' ,'^*{<f')) — C(2), can be studied, also by mean of the exact results 
by Nair, Prabhakar and Sharma PHEll (this being the subject of a paper in preparation 36 ). 

Then, it can be seen that similar results arise from the analysis of a further variant of the 
problem, that we can call "One-in-Four" Matching. We restrict the ensemble of valid matchings to 
permutations n such that, for each i — 1, . . . , n, either 7r(2i — 1) or 7r(2i) are in {2i — 1, 2i}, i.e. in 
the graphical language of section El we have four *s in the four entries of the 2x2 diagonal blocks. 
For this problem, we can describe the configurations by mean of triplets {a, r, n), with (?, r e {0, 1}" 
and TT G Sn with no fixed points, and we can repeat the argument above, such that the best over tt of 
(ij, r) pairs differing also for a single entry are strongly decorrelated, as the sub-instances differ by a 
whole row or a whole column. Nonetheless, Onc-in-Four Matching is trivially reduced to Matching, 
and the same holds for Assignment. Indeed, consider now classes of configurations (ct, t, tt) sharing 
the same tt. The set of 2^" costs for various (a, r) is 



I Yl ^2i~a{i),27r{t)~r(7r{t))\_^ ^^ 



re{0,l}" 

and is trivially minimized, as each variable ai or tj enters only once in the cost function 

--Sn^^-(5I^2^-<T(^),27^(^)-T(7^(^))) =51 mjn «^2»-^,2^(»)-r (31) 

% I 

that is, the minimization problem is reduced to a traditional assignment problem, where 2x2 non- 
diagonal blocks are replaced by the minimum weight among the four entries, and the diagonal blocks 
by infinite weights. 

The same construction can not be repeated for One-in-Two Assignment, as the spin variables 
enter in the cost function in different points. 

Another derivation of this property is given by the reduction of One-in-Four Assignment to 
a Min-Cut-Max-Flow Problem (cfr. 37 for all definitions), which is a polynomial problem. In a 
similar fashion to the mechanism above, the reduction can be done because of the factorization of 

J In the sense of probability, that is, for N the size of the problem, always up to a probability pjv, such that pjv ~^ 

for N -^ oo. 
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all the edges have capacity 1, all the edges except the ones in the middle layer have cost zero, and, 
for the middle layer, an edge connects the i-th vertex on the left to the j'-th vertex on the right only 
if \i/2] 7^ ri/2] , and in this case it has cost Wij . Checking that this instance of Min-Cost-Max-Flow 
reproduces the One-in-Four Assignment, where the maximum flow is n, is a simple exercise, while 
the procedure of equation l|lTT)l corresponds to the "resolution" of the Min-Cost-Max-Flow problem 
on the subgraphs 

min(^ij ) 



6. A Quantum Adiabatic Algorithm for One-in-Tv^ro Matching 

Here we describe a quantum Hamiltonian reproducing a Quantum Adiabatic Algorithm for One-in- 
Two Matching or Assignment, which is a very natural one for condensed-matter systems. We assume 
our instance is composed of a matrix of size 2n, with weights Wij out of the diagonal blocks, all 
much larger than zero (this is always the case, up to a trivial gauge transformation on the instance). 

Consider two layers, each consisting of 2n x 2n discrete positions on a square lattice. We call 
them the row- and the column- layers. We should imagine them geometrically "facing", like two 
plates of a condenser. 

In each row of the row-layer there is one particle, which could jump from one column to a 
neighbouring one. If particle in row i is in layer position (ij), we say that it is in the state |j),./j-). 
Introduce the set of quantum oscillators a^j-, a'l . 

Similarly, in the column-layer we have one particle per column, which can hop on neighbouring 



cU) 



rows, and the state of column j is 
quantum oscillators for this layer. 

A kinetic "diffusion" term would have the form 

Hamiltonian, Hq, is 



if the particle occupies position i. We call 6, 



bl 



the 



a;,i±ia' 



ij- 



So the purely kinetic part of the 



i?n = 



i,3 



[^zJ+l'^hJ 



■j,i-i"«j 



&],,+ !&,.. + 6],,_ifo. 



J,v 



(32) 



On the boundaries of the chains, you can assume either open or periodic boundary conditions (the 
first one is more "physical" , the second one makes the analysis simpler, anyhow the difference is 
negligible). Then, under this Hamiltonian, the particles are non-interacting. At sufficiently low 
temperatures, scaling with n~^, the single-particle system (say, of row «), up to an exponentially 
small fraction of the time, will occupy the ground state with zero momentum, (2n)~^/^ ^ - |j%(j), as 
it has a polynomially-small but finite gap, given by the square of the smallest non-zero momentum 
on the lattice. 

The ground state of the system under the Hamiltonian Hq is 



1^^)03 = (2n)-'"(|l) + --- + |2n)) 



«)(|1) -!-■•• + 



\2n)' 



r{2r> 



2n) 



c(l) 



®(|1) 



\2n)] 



c(2n) 



(33) 



One-in-Two-Matching Problem is NP-complete 20 

and it is a reasonable initial state, for a quantum adiabatic paradigm in which the coordinate-space 
representation reproduces the set of feasible solutions, although it is probably a poor initial state for 
a traditional quantum computing approach, as it is a "tree state" with 0{n^ Inn) symbols (cfr. |2()j . 
ch. 13). 

Up to now, we have one particle per row on the row-layer, and similarly on the column-layer, 
but we have not implemented neither the matching constraint, that row-particle i is in column j iff 
column-particle j is in row i, nor the one-in-two constraint. 

At this aim we introduce an interaction Hamiltonian Hi , which couples the particles on the two 
layers. Furthermore, we will have special terms on the diagonal 2x2 blocks. Introduce the shortcut 

{a''aWh)i^j := a\jai,jb]J)j^i ; {a''aWb)i := (a*a6*6)i,i . (34) 

For pairs {ij) not in a diagonal block, we would have a term of the form —Wij{a'*ab^'b)ij, while 
for a block on entries {2i — l,2i} we would have a combination of the form W (^{a^ab''^b)2i — 
cV(a*a 6*6)21-1 — c), with c a value between and 2 (for example, 1, but notice how the fine- 
tuning of c exactly equal to 1 is not required), and W an energy value larger than the typical 
Wij . Remark that the procedure works also if W and c are slightly varying from block to block, as 
expected from an experimental realization. 

Summarizing, the Hamiltonian Hi (notation: ^^ ■ means "pairs {ij) not in a diagonal block" , 

i.e. such that [|] / [§]) 

n 

Hi = -Y^w.j {a' a b'b\j -WY, {{a'ab'b)2r - c) ((a+a 6*6)3^-1 - c) (35) 

i,j i=l 

is diagonal in coordinate space, as it only involves "number" operators, and has the optimal one-in- 
two matchings as ground states. Indeed, violating a matching constraint would have a finite positive 
cost ~ Wij , while violating a one-in-two constraint would have a cost ~ W. Finally, for a valid 
matching, the energy of the state reproduces exactly the cost function. 

For Wij not scaling with n (which is the proper choice in order to have an extensive energy) , in 
the case of matching we have a gap of order 1, while in the case of random assignment we have a 
gap only algebraically small in n, and not exponentially (cfr. J34II36| ). 

So, the adiabatic interpolation, for time t G [0, 1], given by Ht = {1 — t)Ho + tHi, could be a 
good prescription for a quantum adiabatic algorithm for the one-in-two matching, as, at least at the 
two endpoints, there are no exponentially-small energy gaps. A further analysis of what happens at 
intermediate times would be too long and difficult, and out of the aims of this paper. 

7. Conclusion and perspectives 

In this paper we introduced a problem, the One-in-Two Matching, which lies in between the 
(polynomial) Matching Problem and the (NP-complete) Hamiltonian Circuit Problem (HC). Despite, 
at a first glance, it could seem much more similar to Matching than to HC, we prove that indeed it 
is NP-complete (section |21) . We thus expect that it captures the core reasons of complexity of HC, 
although largely simplifying the model. 

The reduction is relatively simple and specially compact. As a side result, combination with 
a reduction proof from Onc-in-Two Matching to 3-Dimensional Matching (3DM) given in section O 
provides a simple reduction proof for 3DM, which could be preferred for pedagogical purposes to 
the original one by Karp in 1972. 

Everything above applies to the "weighted" variants of the problems: (One-in-Two) Matching 
-^ (One-in-Two) Assignment, Hamiltonian Circuit -^ Travelling Salesman, 3DM — * numerical 3DM. 

In the introduction, and in section we also give very preliminary arguments towards the 
fact that One-in-Two Matching is a viable candidate for an experimental realization of a Quantum 
Adiabatic Algorithm protocol, devoted to the solution of NP-complete problems. In two words, 
the main advantages are the structural simplicity of the problem formulation, and the unusual 
compactness of the reduction. 
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At the aim of still improving this second point for real-life applications (in which short clauses 
could appear, that, when reduced to SAT, would largely increase the instance size), in section^we 
give compact encodings of a larger family of clauses, including SAT, NAE-SAT, One-in-Three-SAT, 
Colouring and much more. Furthermore, in |Appendix A| we find an encoding for all clauses involving 
up to four literals. 

Further directions of investigation will include a better understanding of the statistical 
properties of Random One-in-Two Assignment, at the light of the results of Nair, Prabhakar and 
Sharma OESj, and also with Cavity Method techniques jHE), in relation to the important question 
of determining the geometrical characteristics of the hard instances of TSP |12[ 15^ . 

Appendix A. Truth-table dictionary 

In this Appendix we classify the encoding of all truth tables up to 4 literals. Given k literals ui, 
. . . , Uk, in principle there are 2^^ ^ possible truth tables, corresponding to the set of subsets of 
{True, False} '^ (the satisfying assignments of the literals). More precisely, a clause Ct is identified 
by the set T of literal configurations evaluated to True 

T == {f} C {True, False}''' ; CT(r) = True ^ teT. (A.l) 

Still in section 0] we discuss the "gauge" invariance which allows to consider truth tables inside the 
same class as equivalent, this fact allows to reduce the number of tables to be analyzed, for a fixed 
value of the length k. This number is further reduced by the fact that certain tables are "trivial" 
inside a boolean decision problem, in a sense which is depicted by the following rules: 

Rule 1: The clause is always UNSAT, i.e. T = 0. The instance is determined to be unsatisfiable. 
Rule 2: The clause is always SAT, i.e. T = {True, False}'^. One can reduce to a smaller instance, 
where this clause is removed. 

Rule 3: The clause can be encoded in logical form with no use of a literal (say, Ui). On the truth 
table, this means that, defining the symmetry operation 

i?j(Tl,...,Tfc) = (Tl,...,T,,...,rfc), 

T — Ri{T). The clause is equivalent to a smaller one, and thus the whole instance can be 
reduced in size. 

Rule 4: The clause forces the assignment of a given literal (say, Ui — True, or Ui = False). One can 
simply regard Ui as a synonimous of the fixed boolean variable, and the instance is reduced in 
size (the new instance has one literal less, and all the clauses in which Ui appears are shorter). 

Rule 5: The clause forces the relative assignment of a given pair of literals (say, UiVuj = True, or 
UiVuj — False). One can simply regard Uj as a synonimous of Ui (or of m^), and the instance is 
reduced in size (the new instance has one literal less, and all the clauses in which both Ui and 
Uj appear are shorter). 

Actually, Rule 5 is not so rigorous: sometimes synonima could be useful in certain constructions, and 
for example we made large use of them in our reduction proof from SAT to One-in-Two Matching 
(and we also described the "consistency check" clause, which indeed violates Rule 5 on all the pairs). 
So, in our enumeration we will include all clauses violating Rule 5, although in this case we will 
point out the fact. 

We make a complete enumeration "by hand" up to length 3, where the truth tables are easy to 
handle. We skip both clauses which are trivial w.r.t. rules 1-3, and clauses of the kind "range-T", 
for which a general recipe has already be given in section^ 

For length 4, the manual enumeration would have been a formidable task, so we preferred to 
use a randomized search method: we just tried to fill randomly matrices W, with a given density, 
and classified the resulting clauses, up to when each class had at least one representative. ft 

ft Although this was not assured to happen, we were lucky, and found a representative per each class in a few hours 
of computer time. 
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Another difRculty is that, as the same truth tables are difficult to handle, our results would be 
hardly fruible if the user would not have a recipe to identify his clause of interest in our list. For this 
reason, we decided an easy algorithm for the determination of an "algebraic signature" of the gauge 
orbit of a given truth table. This signature is an integer number, and thus, as in a "dictionary" , 
decides the sorting of the entries of our enumeration. For fruibility to the user, also the clauses 
violating rules 3 and 4 are listed, with a side annotation. On the other side, the cases \T\ < 2 and 
\T\ > 14 = 2^ — 2, which fall in the general framework of section^ have been omitted. 

More detailed explanations are given below. 

Appendix A.l. Clauses of length 2 

The cases \T\ = 0, 1,4 are trivial. The case |r| = 3 is SAT. The case |r| = 2, which is both NAE- 
SAT and XOR, i.e. truth-setting structure, is already in the framework of section El So, there are 
no new cases. 

Appendix A. 2. Clauses of length 3 

The cases |T| = 0, 1, 8 are trivial. The cases \T\ — 2, 6, 7 have been studied in sections |21 and 0] 

Now consider the case \T\ = 3, and call T — {fi,T2,T^3}- The set of "distances", where the 
distance is defined as d{T^ f') = #{i : r^ j^ r/}, can form the triangles (1, 1, 2), (2, 2, 2) and (1, 2, 3). 
The clause (1,1,2) violates Rule 4, because leaves a sub-hypercube (just a square in this case) 
outside. The clause (2, 2, 2) is the One-in- Three clause, which is a special case of range-T clause. 
The clause (1, 2, 3) violates Rule 5. An example is (uiVu2) A (ui Vu2 VM3), where violation of Rule 5 
is made explicit. Anyhow, we give a matrix encoding in Table \Kl\ 

The case \T\ = 5 has a similar case study: if we consider the complementary of T, it has 
cardinality three, and with the same reasoning above we have the three cases with sets of distances 
(1, 1, 2), (2, 2, 2) and (1, 2, 3). The difference is that now all these cases are new. 

In the case \T\ = 4, either all the four vectors are at distance 2 (and in this case we deal 
with the important XOR clause) or there are at least two neighbouring vectors (say, ti^2)- Then, 
six vertices remain on the cube, of which four are neighbours of ti^2, so either -73^4 fill the pair of 
non-neighbouring sites, and we have the trivial clause (M1VU2), which violates many rules, or there 
is a vertex with two neighbours, say, ti is at distance one from vz and 73. If it is at distance one also 
from r4, we have an "At-least-2-in-3" clause, i.e. a special case of range-T clause with k = hmax = 3 
and /imin = 2, for which we have a general recipe. Otherwise, if T4 is at the last vertex of the 
square with vertices f 1^2.3, we have a trivial clause of length 1. There are only two other choices 
left, either the pairs of neighbouring vertices form an open path of length 3 (the "first-F-then-T" 
clause, (ui A U2)V(u2 A U3), with truth table {{TTT), (FTT), (FFT), (FFF)}), or the last vertex 
T4 is antipodal to ri . 

This completes the case study for fc = 3: we have three new cases at |r| — 4 and three at 
|r| — 5, and a case at |T| = 3 which violates Rule 5. The suggested encoding matrices are listed in 
Table EH 

Appendix A. 3. Clauses of length 4 

Each entry of our dictionary is encoded in the form 

|r| Algebraic signature Truth table Encoding matrix W #{wij = 1} (Comments) 

The Algebraic signature allows to determine the gauge class of a given truth table. This is the reason 
why the truth tables are sorted according to that. In a sense, it is a univoque name for a clause: if 
you want to encode a given truth table with 4 literals and |r| True assignments into a One-in-Two 
pattern, and this truth table does not follow into one of the general classes studied in sections El 
and01 you should calculate the algebraic signature (a simple and fast algebraic procedure), and then 
search the corresponding entry in this dictionary; then, the other fields of the entry will allow you 
to construct the proper matrix. 
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\T\ 



Table Al. Encoding matrices for the new clauses of length 3. For each clause are given: a 
boolean description, a suggested One-in-Two encoding, and a graphical description of the truth 
table (bottom vertex: (TTT); bottom row, left to right: (TTF) ■ ■ ■ (FTT); filled disk: r G T.) 
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Finding a good algebraic signature recipe is in principle a complicated task, because asking 
whether two truth tables are isomorphic subgraphs of the fc-dimensional hypercube is indeed a 
graph- isomorphism problem, which is a hard problem (cfr. jSH])- Nonetheless, in our case it is 
certainly affordable, because the coordination is bounded by 4. Furthermore, as fortunately the 
hard examples of graph-isomorphism are only a narrow set with special geometrical properties, it is 
plausible that also the naive algebraic techniques, consisting of a calculation of a determinant, allow 
to discriminate the classes. Consider the distance d{T,T') = #{z : Ti ^ r,'}, which is an integer in 
{0,1,..., k}, and only for coincident vectors; then consider fc -I- 1 relatively-trascendent real values 
r{i). We have chosen 



r(0) 
r(l) 
r(2) 
r(3) 
r(4) 



= 1; 

= C(3)"^ = 0.831907372580707469; 
= -7Euicr = -0.577215664901532861 ; 
^ e-i = 0.367879441171442322 ; 
= -TT-i = -0.318309886183790672. 



(A.2) 



The choice is, of course, arbitrary. We choose decreasing real numbers in [0,1], with alternating 
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signs, in order to have determinants of the same order of magnitude. Then, we build the square 
matrix M of size \T\, with Mf^f^/ = r(d{f,f')). Just to reduce computation times, if \T\ > 2*^^^ we 
use instead of T its complementary set. The chosen "algebraic signature" of truth table T is the 
integer part of 10^^ det M{T). 

The truth table is simply encoded in a compact way, the 2^^ possible vectors r being ordered 
lexicografically, and the table T being denoted by the integer n{T) G {0, . . . , 2^^ — 1} corresponding 
to the binary string of vectors evaluated to True 



n(r) - J2 2 



n'{T) 



tGT 



n'(T)=^2'^-^5,.,TVue. 



(A.3) 



The representative of the class is chosen in arbitrary way. Indeed, it comes out that it is the one 
with lowest index niT), among the Ts in the same class. 

Similarly, the encoding matrix is written in a more compact way as a vector of integers, where 
each row of Os and Is, of length 2fc = 8, is replaced by the corresponding integer in {0, . . . , 2® — 1}. 
Just for statistic reasons, the number of Is in the matrix is indicated. 

Finally, if the clause is trivial in the sense of the rules 3-5 above (failures of rules 1 and 2 are 
excluded trivially), we report the fact, and the rule it fails (if it is trivial for more than one reason, 
we choose a reason of smaller rule index) . 

Here the dictionary follows: 
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U3 <^C 


4 


-55489630922495513 


8242 


(3,1,128,192,24,0,8,48 


11 


R4: 


Us — True 


4 


-45874752110564406 


8227 


(8,8,3,16,161,0,128,96 


11 


R4: 


Us — True 


4 


-37719695731007326 


8211 


(40,8,136,16,3,0,128,96 


11 


R4: 


Us — True 


4 


-32798143896572570 


8203 


(9,8,128,64,3,4,32,48 


11 


R5: 


ui\/U2 — False 


4 


-30712541418824951 


8199 


(32,16, 132,1,9,128,2,72 


11 


R5: 


ui\/U2 — False 


4 


-28706142641542548 


8338 


(3,8,32,192,24,16,8,133 


13 






4 


-28647676317807325 


8481 


(8,64,128,8,48,0,2.1 


8 


R3: 


Ui^C 


4 


-26253185098326629 


8205 


(6,16,34,64,8,128,160,1 


11 


R5: 


ui\/U2 — False 


4 


-23870314469588248 


8218 


(3,64,12,16,129,32,32,136 


12 






4 


-21717939174696339 


8214 


(10,64,129,8,144,4,32,33 


12 






4 


-21459144106706130 


8229 


(6,64,128,16,9,32,34,8 


11 


R5: 


ii2Vii4 — False 


4 


-19374761209641485 


8217 


(33,64,160,1,8,128,2,20 


11 






4 


-17236351616013344 


8220 


(128,16,5,64,9,32,2,40 


11 


R5: 


ii2Vii3 — True 


4 


-15085506227251376 


8580 


(32,64,128,16,1,4,2,8 


8 


R5: 


U2'^U4 — False 


4 


-14842292216336441 


8772 


(2, 1,128,16,8,4,32,64 


8 


R3: 


U3^C 


4 


-13941487138238279 


8265 


(12,1, 129,16,96,128,2,40 


12 






4 


-9019935303803522 


8520 


(4,24, 160,129,64,17,2,40 


13 


R5: 


uiVu^ — True 



5 -137319486169585905 8370 

5 -131925251824472796 8243 

5 -106990011212970741 8235 

5 -106200686484487583 8467 

5 -100711446287411196 8483 

5 -97383480909282057 8207 

5 -93806950451816800 8215 

5 -87252051588588716 8231 

5 -84213109538527750 8219 

5 -81265445872985527 8339 

5 -79740217163681572 8269 

5 -72810579632836741 8221 

5 -69543625692646563 8222 

5 -69010625791285801 8342 

5 -67084531825011153 8267 

5 -65196296191649679 8252 

5 -64617744063353716 8237 

5 -62513764335079739 8297 

5 -60805966899451609 8282 

5 -58556655299675617 8581 

5 -57613278077157236 8582 

5 -54873765588920526 8524 

5 -53148596289303402 8773 

5 -51297235131455269 8476 



(7, 64, 32, 9, 144, 128, 16, 24 

(3, 16,32,17,8,0, 128, 192 

(9, 1, 128, 192, 10,32,32,52 

(9,64,32, 129,24,0,2, 136 

(34,64,160,16,129,0,8,8 

(32,8, 128, 1,24,4,2,65 

(11, 1, 128, 144, 40, 32, 22, 72 

(56, 64, 132, 144, 33, 128, 2, 40 

(130, 16, 9, 72, 49, 4, 128, 40 

(24,1, 129,28,168,1,2, 192 

(34,1, 128,80,8,161,6,20 

(34, 1, 128, 144, 40, 3, 28, 68 

(128, 16, 11,8, 17,4,34,96 

(6, 10, 4, 16, 129, 132, 104, 65 

(160, 1,9,16,96,4,2, 136 

(36,64,128,16,8, 1,2,33 

(161,64, 12, 16,10, 130, 128,33 

(37, 1,6, 18,8,144,128,224 

(34,64,36, 16,8,128, 128, 1 

(128,8,32, 1,66,4, 10,16 

(34,1, 132, 16,8,32, 128,65 

(4,8,3,66,48,128,80,33 

(2,64, 160, 16,8, 128,36, 1 

(3,64,32, 129,16,20,6, 136 



13 
10 
13 
11 
10 
10 
15 
14 
14 
15 
14 
15 
13 
15 
12 
10 
14 
15 
10 
10 
11 
13 
10 
13 



R4: Us = True 

R4: Us = True 
R4: Us = True 
R5: U1VU2 — Fals 



R5: M2VM3 = Tru 



R5: U2VU4 — False 

R5: uiVus — True 
R5: ui\/U4 — False 



One-in- Two- Matching Problem is NP-complete 
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5 


-49945641870770219 


9289 




5 


-49144516830228831 


8492 




5 


-48332850287105365 


8521 




6 


-264953748043115284 


8755 




6 


-255903677087804212 


8371 




6 


-246423116736164502 


8499 




6 


-234055909760095663 


8251 




6 


-221848639666228098 


8247 




6 


-219506781223121329 


8374 




6 


-211897911062832022 


9009 




6 


-198829802885264160 


8239 




6 


-196931604634552768 


8223 




6 


-193105575805408406 


8314 




6 


-186798665756457220 


8271 




6 


-186227373079063234 


9293 




6 


-184055504850097333 


8299 




6 


-180868489214616524 


8254 




6 


-176493283894497754 


8343 




6 


-174574944498457624 


8347 




6 


-164368307347587140 


8477 




6 


-164359755773378736 


8283 




6 


-159842111719858328 


8253 




6 


-157236364066051777 


8583 




6 


-155304537441776442 


8286 




6 


-153275171299567776 


8525 




6 


-153237592690580239 


8494 




6 


-150450270880229517 


8301 




6 


-147186103159461859 


8652 




6 


-144210284492927278 


8598 




6 


-142969789437608309 


8478 




6 


-142370542521970898 


8523 




6 


-141698008727709048 


8789 




6 


-138524570147284945 


8749 




6 


-138108485317911594 


8775 


(1 


6 


-135427674056010516 


8589 




6 


-133261805405082912 


8493 




6 


-132786334798600321 


8526 




6 


-130731773807879766 


8781 




6 


-130024582675614020 


9291 




6 


-125686589014661632 


10644 




6 


-125660450966743013 


8857 




6 


-124616208289630152 


8604 




6 


-123666407520778105 


8508 




6 


-122864364540640491 


8613 




6 


-120544022323974311 


8538 




6 


-119127428200296783 


8790 




6 


-118936297852452994 


9321 




6 


-113571267670065141 


8796 




6 


-112145695798823883 


8556 




6 


-110328833598895423 


9306 




6 


-109077852951287423 


9036 




6 


-101641038122898375 


9576 




6 


-101300461292717132 


9156 




7 


-489610810524822682 


9011 




7 


-452337729517226882 


8379 




7 


-408354319245236969 


8375 




7 


-400582064369183280 


8315 




7 


-393459705592617740 


8255 




7 


-384942911261086050 


8751 




7 


-378955175228864660 


8479 




7 


-372730902319860600 


8287 




7 


-365307444265761614 


8351 




7 


-350325235410304638 


8599 




7 


-346991515368820537 


9295 




7 


-345602958711606536 


8495 




7 


-342028055323505620 


8303 




7 


-334343815080691876 


8411 




7 


-331245387788784281 


8527 




7 


-329455080315344929 


8591 




7 


-326446789752148616 


8318 




7 


-309905349059268324 


8539 




7 


-305364326870046401 


8510 




7 


-303440793052772106 


8783 




7 


-303164240726434633 


8615 




7 


-294782498316163996 


8509 




7 


-293737612421108058 


8654 




7 


-291593988308626710 


8606 




7 


-286206440273200448 


8791 




7 


-282360887598956176 


9323 





(136, 33, 131, 10, 96, 4, 82, 48 
(128,40,14,9,18,32,68,65 
(48,40, 136,66,3,4, 192, 17 

(24, 16,32,8, 192,0,2, 1 

(128, 16,1,24,72, 1,2,37 

(29, 12,8,81,65, 133,2, 144 

(40,8, 136, 16, 161, 1,2, 196 

(40, 64, 10, 24, 3, 128, 132, 144 

(6, 64, 128, 16, 129, 130, 44, 9 

(8,8,145,16,103,0,228,96 

(33, 1, 128,8,24,4,2,225 

(11, 16, 136,8, 129,4,48,96 

(3, 1, 160, 16,96,128,8, 12 

(38,1,128,16,8,3,164,68 

(6,5,131,66,8,160,112,48 

(14, 8, 32, 192, 18, 17, 132, 129 

(45, 64, 128, 16, 129, 32, 2, 12 

(12, 10,2, 16,66, 164, 192, 1 

(160, 1,138,16,98,1,8,13 

(130, 1,9, 16,96,4, 128,40 

(129, 1, 36, 74, 146, 132, 24, 40 

(35,1, 128, 16, 194,32, 12, 12 

(26, 16, 160, 72, 161, 4, 138, 40 

(142, 8, 128, 9, 66, 36, 48, 17 

(128,10,9,66,65,36,48,20 

(162, 66, 12, 8, 65, 129, 32, 144 

(49, 64, 136, 129, 18, 34, 132, 12 

(4,8,3,3,80,48,160, 192 

(18,68,4, 11,16,224, 144,5 

(38, 80, 2, 73, 24, 130, 132, 16 

(3,64, 136,129,48,4, 18,40 

(32,8, 128, 16,72,4,2,1 

(37,33,14,8,88,32,128, 145 

(131, 16, 132, 136, 72, 32, 50, 129 

(48, 64, 131, 17, 9, 2, 132, 136 

(128,48,69, 1,8,32,2,84 

(20,48, 130,68,8,34, 192,1 

(148, 64, 33, 136, 16, 130, 4, 33 

(136, 40, 34, 66, 48, 68, 148, 1 

(10, 96, 96, 10, 17, 20, 132, 129 

(6, 1,32, 16,8,128, 128,65 

(128, 18,4,26,9,48,68,65 

(20,64,3, 144,1,34,4, 136 

(40, 16, 128,64,40,2,6,1 

(145, 64, 32, 10, 3, 132, 12, 144 

(128,16,12,8,67,33,18,33 

(36,3,13,66,10,208,160,20 

(33, 18, 36, 192, 72, 128, 10, 17 

(6, 64, 130, 17, 40, 129, 12, 48 

(17,18, 160,192,65,4, 10,40 

(20,8,65,5,96,128,2,48 

(20,80,32,10,3,132,136,65 

(6, 48, 131, 17, 64, 160, 12, 24 

(8,64,32,144,25,0,2, 129 

(128,8,32,82,24,4, 16, 1 

(20, 16, 32, 13, 136, 192, 2, 68 

(11,1, 128, 192,48,4, 18,41 

(136, 1,32,192,24,4,2,9 

(33, 1,136,8,192,4,2,48 

(144, 1,32,9,193,4,2,24 

(10,64, 160, 16,129,4,32,9 

(128, 12, 32, 11, 16, 20, 66, 80 

(6,68,4,24,24, 162, 192,1 

(36,3,192,16,10,34,140,65 

(21, 1, 32, 8, 194, 128, 130, 28 

(168,8,33, 17,144,4,2, 192 

(136,64,32,9, 131, 1,8,21 

(13,72,36, 192,1,52,2, 144 

(129,8,40,1,96,4,2,144 

(130,8,134,16,11,36,40,65 

(150,64,9,3,152,48,36,136 

(19, 64, 136, 10, 176, 36, 128, 5 

(11,64,32,9,146,4, 128,24 

(40, 16, 168, 192, 137, 1, 2, 21 

(144, 8, 36, 129, 72, 33, 2, 48 

(3,24,96,68,40,144, 130,5 

(3,68,4,146,8,48,208,5 

(10,64, 160, 1,40,4, 128, 17 

(129,40,140,16,33,6,70,96 



17 
15 
15 

9 
12 
18 
14 
14 
14 
17 
13 
14 
12 
14 
17 
15 
13 
14 
15 
12 
17 
14 
17 
16 
15 
15 
16 
14 
16 
16 
14 

9 
17 
16 
15 
13 
14 
14 
16 
16 
10 
15 
13 
11 
15 
14 
18 
15 
15 
15 
13 
15 
16 

11 
11 
14 
15 
12 
12 
13 
12 
14 
15 
16 
15 
14 
14 
16 
12 
16 
18 
16 
14 
16 
14 
16 
16 
12 
17 



R3: «3 ^ C 
R4: Us = True 

R4: M3 = True 



R5: uiVus — True 



R3: «3 ^ C 



R3: «3 ^ C 



R5: U2VU4 — False 



R5: uiVus — True 
R5: U1VU3 — True 
R4: M3 = True 



One-in- Two- Matching Problem is NP-complete 
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7 


-281531220380667350 


8653 




7 


-280171463154608827 


9325 




7 


-280099211892526104 


10646 




7 


-276166021834513026 


8765 




7 


-276136614689952116 


9310 




7 


-275318071137808221 


8605 




7 


-270984298897396969 


9307 




7 


-268713156635853130 


8542 




7 


-263998281304105199 


8797 




7 


-261677808517195389 


10645 




7 


-259011690208036504 


8558 




7 


-258182304328415778 


8557 




7 


-251249744669286557 


9037 




7 


-249858540374917035 


8861 




7 


-247427194046103816 


9622 




7 


-247143106478950578 


8621 




7 


-243566351808186474 


8813 




7 


-242068486226572112 


8798 




7 


-237212856195456665 


9628 




7 


-227914280620238644 


9577 




7 


-226039824958208568 


9580 




7 


-223022686093957866 


9578 




7 


-221577833850858925 


9116 




7 


-217580664595077813 


9562 




7 


-217251409270906367 


9673 




7 


-213455902944141485 


9158 




7 


-211010290780943168 


9157 




7 


-206712752788130826 


9052 




7 


-205559919351541300 


9164 




8 


-899283463066040377 


13107 




8 


-772207269962157112 


8443 




8 


-721019390851251273 


8383 




8 


-686263475047731580 


8415 




8 


-673241419621910315 


8319 




8 


-666032231320527476 


8767 




8 


-661718979480894926 


8511 




8 


-650672248853117917 


8631 




8 


-642756855300046707 


8571 


(1 


8 


-636308968528110613 


8543 




8 


-613284019874591394 


9311 




8 


-602996490440398480 


8607 




8 


-583125401929517341 


8863 




8 


-578102585724439884 


8623 




8 


-566727496236131365 


9339 




8 


-564929920852814672 


8823 




8 


-560433615140304576 


9623 




8 


-558283470099177545 


10647 




8 


-553542914918976610 


8559 




8 


-553440880392537039 


8655 




8 


-550329284259922678 


8815 




8 


-540872899944542692 


8670 




8 


-537968997943824871 


9039 




8 


-525172959641607591 


8799 




8 


-521929389664079637 


9677 




8 


-521206989824199314 


9327 




8 


-520079249630644275 


27030 




8 


-513543929637538094 


8574 




8 


-510932219079709161 


9435 


(i 


8 


-508220186566127504 


8893 




8 


-503779283191348472 


9118 




8 


-494343963011923576 


8830 




8 


-493714312721210462 


8637 


(5 


8 


-489334924266067936 


8925 




8 


-489101200908942974 


9629 




8 


-488708077465335931 


10709 




8 


-482042543996118538 


8829 




8 


-480275171830561668 


9581 




8 


-479323869615069659 


10710 




8 


-477519058171793838 


9566 




8 


-473957427851750723 


8685 




8 


-471341320523760961 


9579 




8 


-464534444202425179 


9630 




8 


-463299397530769003 


9053 




8 


-457924647789240028 


10653 




8 


-452601408428460382 


9342 




8 


-442700287585267532 


9117 




8 


-441281571114611767 


9706 




8 


-431649046291600881 


9594 




8 


-428443387129934507 


9159 




8 


-425985466039421809 


9133 





(48,72,6,20,40, 130, 192, 1 

(6, 17, 32, 208, 72, 131, 20, 24 

(38, 3, 40, 136, 81, 20, 136, 193 

(132, 1, 160,9,24,33,2,84 

(132,72,32,3,18,5,80,136 

(10,64,6,16,136,32, 160,1 

(132,5,40, 136,96,81,2,20 

(11, 33, 130, 152, 16, 36, 192, 69 

(50,64,4,17,8,130,36,129 

(40, 64, 3, 10, 144, 132, 20, 161 

(5, 20, 160, 129, 24, 3, 80, 104 

(52, 16, 38, 192, 8, 130, 132, 1 

(5,72,36, 17,18,160,10,192 

(160,8,10,16,96,4, 130,1 

(25, 5, 130, 193, 136, 98, 36, 20 

(5,66,32,9,136,128,66,17 

(6, 24, 32, 192, 16, 131, 14, 17 

(5,80, 130,192,56,32,10,9 

(20, 48, 5, 129, 72, 98, 18, 136 

(17,96,6,11,136,132,68,48 

(129, 10, 68, 16, 72, 34, 36, 129 

(20,64,3,144,33,36,10,136 

(20,80,32, 192,8,3,4, 129 

(128,72,96, 10, 18,20,36, 1 

(160, 40, 130, 18, 192, 84, 12, 1 

(20, 64, 160, 140, 1, 34, 132, 25 

(6, 72, 7, 144, 33, 161, 196, 24 

(20, 20, 32, 129, 8, 67, 18, 192 

(10,64, 130, 16, 136,4,32, 1 

(8,64,32,16, 136,0,2, 1 

(128,64,9,16,3,4,32,40 

(130, 64, 32, 13, 129, 6, 144, 16 

(130, 80, 128, 11, 17, 4, 98, 16 

(3, 8, 12, 16, 200, 36, 160, 129 

(132, 1, 160, 16,8,32,2,65 

(132,5,32,81,8, 144,2,144 

(128,17,32, 13,16,84,2,68 

(128, 237, 38, 120, 80, 67, 152, 33 

(130,64,9,16,33,4,128,40 

(128,33,65,68,10,6,24,48 

(9,64,32,131,24, 17,6,20 

(42, 16, 9, 17, 136, 4, 130, 208 

(128,8,34,16,65,4,32,9 

(11, 12, 160, 192, 33, 68, 18, 48 

(48,66, 176, 132,8, 128,2, 1 

(18, 5, 32, 138, 72, 69, 144, 20 

(7, 80, 32, 13, 136, 144, 18, 192 

(137, 64, 128, 24, 131, 4, 34, 48 

(136,40,68,17,192,5,2,48 

(5,64,32, 129, 16, 148,2, 12 

(128,9,69,80,56,3,34,12 

(48, 18, 192, 129,65,4,34,9 

(128,8,34, 16,72,2,36, 1 

(36,65,6,16,8,130,160,65 

(144, 24, 32, 92, 66, 130, 20, 1 

(41, 22, 193, 194, 161, 82, 28, 44 

(144, 20, 65, 5, 136, 164, 2, 104 

(89, 42, 150, 66, 114, 194, 14, 193 

(22, 20, 65, 208, 8, 162, 36, 129 

(128, 48, 34, 67, 24, 48, 76, 5 

(35, 64, 44, 144, 40, 129, 20, 129 

(28, 64, 160, 129, 138, 33, 10, 144 

(130,8,32,16,66,4, 136,1 

(130, 48, 96, 12, 66, 20, 136, 1 

(42, 64, 12, 20, 138, 34, 144, 1 

(3, 3, 128, 16, 136, 68, 44, 192 

(6,5,32,208,8,130,20,65 

(57, 6, 201, 71, 161, 82, 92, 172 

(48,65,5,192,9,36,2,152 

(138, 64, 34, 137, 48, 4, 144, 129 

(11, 96, 36, 12, 3, 130, 208, 48 

(17, 48, 192, 137, 129, 6, 6, 104 

(132,20,32,82,1,130,80,9 

(41, 64, 129, 18, 129, 4, 18, 40 

(51, 64, 10, 144, 161, 4, 160, 136 

(36,64, 129, 16,8, 128,2,33 

(5,64,36,129,8,50,6,144 

(128,72,3,10,48,4,80,33 

(144, 20, 32, 136, 65, 6, 66, 17 

(176, 33, 13, 144, 10, 160, 34, 69 



15 
17 
19 
15 
15 
12 
16 
18 
14 
16 
17 
15 
16 
12 
19 
14 
16 
16 
17 
17 
15 
15 
13 
14 
16 
16 
19 
15 
11 

8 
11 
14 
15 
15 
11 
14 
14 
24 
12 
15 
15 
16 
11 
17 
13 
17 
17 
15 
15 
13 
17 
15 
11 
14 
16 
24 
17 
26 
18 
17 
17 
17 
11 
15 
16 
15 
15 
28 
16 
16 
18 
18 
15 
15 
17 
11 
15 
14 
15 
19 



R5: uiVus — True 
R3: «i ^ C 



R3: «3 ^ C 



R3: «3 ^ C 
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8 -424439857990218041 9675 

8 -421336022871598134 9654 

8 -417702864512524473 9054 

8 -413621507243457915 9708 

8 -410176819348621982 9166 

8 -409845545228962591 9690 

8 -405077671322438068 9582 

8 -398409201621308822 10104 

8 -393313530918741654 11730 

8 -393246810197863372 9174 

8 -387574329040405188 9084 

8 -382828422517523851 9165 

8 -380718107299096901 9705 

8 -376733361301144475 15555 

8 -370347967947654772 9180 

8 -369336898403981540 9189 

8 -369049161187396959 9660 

8 -347268859906627992 11220 

8 -344304718756045964 10200 

8 -318880083885350459 10212 

8 -284276962660566950 13260 



9 


-489610810524822682 


56524 




9 


-452337729517226882 


57156 




9 


-408354319245236969 


57160 




9 


-400582064369183280 


57220 




9 


-393459705592617740 


57280 




9 


-384942911261086050 


56784 




9 


-378955175228864660 


57056 




9 


-372730902319860600 


57248 




9 


-365307444265761614 


57184 




9 


-350325235410304638 


56936 




9 


-346991515368820537 


56240 




9 


-345602958711606536 


57040 




9 


-342028055323505620 


57232 




9 


-334343815080691876 


57124 




9 


-331245387788784281 


57008 




9 


-329455080315344929 


56944 




9 


-326446789752148616 


57217 




9 


-309905349059268324 


56996 




9 


-305364326870046401 


57025 




9 


-303440793052772106 


56752 




9 


-303164240726434633 


56920 




9 


-294782498316163996 


57026 




9 


-293737612421108058 


56881 




9 


-291593988308626710 


56929 




9 


-286206440273200448 


56744 




9 


-282360887598956176 


56212 




9 


-281531220380667350 


56882 




9 


-280171463154608827 


56210 




9 


-280099211892526104 


54889 




9 


-276166021834513026 


56770 




9 


-276136614689952116 


56225 




9 


-275318071137808221 


56930 




9 


-270984298897396969 


56228 




9 


-268713156635853130 


56993 




9 


-263998281304105199 


56738 




9 


-261677808517195389 


54890 




9 


-259011690208036504 


56977 




9 


-258182304328415778 


56978 




9 


-251249744669286557 


56498 




9 


-249858540374917035 


56674 




9 


-247427194046103816 


55913 




9 


-247143106478950578 


56914 




9 


-243566351808186474 


56722 




9 


-242068486226572112 


56737 




9 


-237212856195456665 


55907 




9 


-227914280620238644 


55958 




9 


-226039824958208568 


55955 




9 


-223022686093957866 


55957 




9 


-221577833850858925 


56419 




9 


-217580664595077813 


55973 




9 


-217251409270906367 


55862 


(5 


9 


-213455902944141485 


56377 




9 


-211010290780943168 


56378 




9 


-206712752788130826 


56483 




9 


-205559919351541300 


56371 





10 -264953748043115284 56780 
10 -255903677087804212 57164 
10 -246423116736164502 57036 



(25, 80, 34, 129, 136, 162, 6, 68 

(12, 72, 130, 80, 129, 38, 160, 17 

(128,3,65,25,16,36,34,76 

(130, 18, 68, 24, 72, 164, 36, 1 

(36,10,136,144,96,20,66,1 

(25,10,132,20,96,160,66,5 

(129, 24, 130, 80, 65, 36, 34, 12 

(160, 48, 68, 17, 10, 130, 12, 65 

(162, 12, 36, 3, 200, 81, 14, 144 

(37, 64, 128, 19, 152, 48, 130, 12 

(17, 64, 40, 18, 33, 132, 130, 12 

(20, 20, 192, 129, 8, 34, 6, 65 

(12, 64, 34, 129, 136, 49, 6, 144 

(35, 80, 67, 224, 92, 102, 195, 112 

(36,64, 128,16,3,6, 168,9 

(133,10,160,72,24,6,96,33 

(19,40,68, 17,9, 132,6,224 

(160, 76, 132, 66, 24, 34, 20, 1 

(138, 24, 69, 14, 67, 193, 100, 60 

(6, 73, 136, 131, 96, 20, 162, 81 

(32,64, 130, 16,8,4, 130, 1 

(3, 1, 160, 16,40,4,8, 192 

(12, 72, 128, 16, 129, 164, 2, 65 

(160,8, 130, 16,74,4,40,1 

(129, 8, 33, 21, 72, 34, 18, 148 

(50, 64, 128, 144, 136, 4, 40, 1 

(12, 64, 34, 1, 136, 128, 32, 17 

(42, 64, 136, 144, 56, 4, 128, 1 

(41,64, 128, 131,24,4, 18, 16 

(21,66,2,193,8,144,38,20 

(19, 80, 36, 129, 8, 129, 34, 196 

(22, 64, 32, 9, 152, 5, 136, 129 

(130, 64, 9, 136, 48, 4, 144, 33 

(144,96,65, 136,9,6,4,49 

(6,3,32,80,136,192, 12,17 

(12, 48, 96, 11, 17, 132, 18, 193 

(20, 80, 32, 131, 8, 18, 134, 192 

(131, 64, 34, 24, 9, 160, 132, 20 

(42, 9, 134, 16, 73, 160, 36, 193 

(136, 72, 32, 133, 18, 17, 66, 68 

(5,18,32,88,17,6,192, 136 

(21,80,32,7,8,82, 132, 144 

(17,12,32, 192,9,68,2, 145 

(6, 35, 68, 144, 24, 192, 14, 65 

(129, 13, 3, 84, 40, 208, 66, 48 

(48, 64, 130, 12, 144, 132, 34, 1 

(3, 12, 191, 88, 210, 41, 83, 186 

(42, 64, 4, 129, 136, 164, 18, 17 

(17,66,32,152,17,4,66,136 

(41, 22, 201, 194, 161, 82, 30, 44 

(132, 80, 34, 66, 40, 20, 136, 1 

(48,5,65,144,9,38,6,200 

(17, 72, 160, 13, 80, 132, 2, 176 

(176, 68, 3, 21, 8, 162, 192, 21 

(25, 68, 36, 13, 129, 98, 130, 144 

(17,64,4, 144,8,34,6, 145 

(21, 18, 168, 194, 161, 6, 104, 24 

(5, 68, 6, 26, 24, 112, 160, 129 

(11,66,32, 136, 17,20,68, 145 

(132, 33, 96, 10, 24, 129, 18, 68 

(129, 64, 33, 136, 49, 4, 2, 24 

(160, 42, 12, 82, 17, 70, 200, 129 

(160, 64, 13, 24, 129, 36, 2, 144 

(17, 40, 192, 192, 24, 34, 6, 5 

(36, 80, 130, 10, 72, 132, 48, 1 

(9,74,96,194,161,5, 18,20 

(5, 18,32,137,24,6,80, 192 

(5,18,36,17,8,34,192,192 

(136,33,65,84,9,52,2,148 

(13, 65, 37, 136, 144, 48, 2, 192 

(130, 18, 96, 137, 24, 4, 48, 65 

(25, 40, 134, 18, 66, 100, 164, 129 

(20, 64, 32, 140, 17, 130, 6, 9 

(5,64,36,136,9,146,34,17 

(20, 33, 34, 144, 72, 5, 66, 137 

(128, 1,33, 16,40,4,2,72 

(128,64,33, 16,8,32,2,5 

(3,64,32,28, 16,6, 144, 129 

(6, 10,32, 17, 192, 132,76, 16 



R3: «i ^ C 



R3: Ml ^ C 



18 
17 
16 
16 
15 
17 
16 
16 
20 
17 
15 
15 
16 
26 
14 
17 
18 
16 
24 
20 
10 

12 
14 
13 
17 
13 
12 
14 
14 
17 
17 
16 
14 
16 
15 
18 
16 
16 
19 
16 
16 
17 
15 
18 
19 
14 
30 
16 
15 
26 
15 
18 
17 
19 
19 
14 
21 
18 
17 
16 
14 
20 
15 
16 
15 
19 
16 
15 
18 
17 
16 
20 
15 
16 
17 
11 



10 R3: «3 

14 

15 
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10 


-234055909760095663 


57284 




10 


-221848639666228098 


57288 




10 


-219506781223121329 


57161 




10 


-211897911062832022 


56526 




10 


-198829802885264160 


57296 




10 


-196931604634552768 


57312 




10 


-193105575805408406 


57221 




10 


-186798665756457220 


57264 




10 


-186227373079063234 


56242 




10 


-184055504850097333 


57236 




10 


-180868489214616524 


57281 




10 


-176493283894497754 


57192 




10 


-174574944498457624 


57188 




10 


-164368307347587140 


57058 




10 


-164359755773378736 


57252 




10 


-159842111719858328 


57282 




10 


-157236364066051777 


56952 




10 


-155304537441776442 


57249 




10 


-153275171299567776 


57010 




10 


-153237592690580239 


57041 




10 


-150450270880229517 


57234 




10 


-147186103159461859 


56883 




10 


-144210284492927278 


56937 




10 


-142969789437608309 


57057 


(1 


10 


-142370542521970898 


57012 




10 


-141698008727709048 


56746 




10 


-138524570147284945 


56786 




10 


-138108485317911594 


56760 




10 


-135427674056010516 


56946 




10 


-133261805405082912 


57042 




10 


-132786334798600321 


57009 




10 


-130731773807879766 


56754 




10 


-130024582675614020 


56244 


(2C 


10 


-125686589014661632 


54891 




10 


-125660450966743013 


56678 




10 


-124616208289630152 


56931 




10 


-123666407520778105 


57027 




10 


-122864364540640491 


56922 




10 


-120544022323974311 


56997 




10 


-119127428200296783 


56745 




10 


-118936297852452994 


56214 




10 


-113571267670065141 


56739 




10 


-112145695798823883 


56979 




10 


-110328833598895423 


56229 




10 


-109077852951287423 


56499 




10 


-101641038122898375 


55959 




10 


-101300461292717132 


56379 




11 


-137319486169585905 


57165 




11 


-131925251824472796 


57292 




11 


-106990011212970741 


57300 




11 


-106200686484487583 


57068 




11 


-100711446287411196 


57052 




11 


-97383480909282057 


57328 




11 


-93806950451816800 


57320 




11 


-87252051588588716 


57304 




11 


-84213109538527750 


57316 




11 


-81265445872985527 


57196 




11 


-79740217163681572 


57266 




11 


-72810579632836741 


57314 




11 


-69543625692646563 


57313 




11 


-69010625791285801 


57193 




11 


-67084531825011153 


57268 




11 


-65196296191649679 


57283 




11 


-64617744063353716 


57298 




11 


-62513764335079739 


57238 




11 


-60805966899451609 


57253 




11 


-58556655299675617 


56954 




11 


-57613278077157236 


56953 




11 


-54873765588920526 


57011 




11 


-53148596289303402 


56762 




11 


-51297235131455269 


57059 




11 


-49945641870770219 


56246 




11 


-49144516830228831 


57043 




11 


-48332850287105365 


57014 


(15 


12 


-64417460349479835 


56797 




12 


-55489630922495513 


57293 




12 


-45874752110564406 


57308 




12 


-37719695731007326 


57324 




12 


-32798143896572570 


57332 





(137, 64, 32, 140, 3, 6, 144, 144 

(50,64,132,144,8,33, 160,5 

(128,35,6,80,8,49,70,68 

(3, 65, 34, 20, 40, 132, 24, 192 

(128, 16,33,80,41,4,2,40 

(18, 64, 32, 138, 144, 148, 6, 1 

(3, 64, 130, 144, 8, 36, 20, 49 

(12, 8, 96, 144, 19, 130, 160, 65 

(136, 7, 133, 74, 48, 50, 164, 65 

(56, 64, 3, 152, 129, 4, 34, 145 

(128, 17,32,22,8,82,20,65 

(3,88,33,140,81,68,130,48 

(6,18,160,192,40,5,88,33 

(12, 72, 33, 17, 144, 38, 130, 208 

(20, 68, 32, 7, 24, 146, 192, 9 

(137, 64, 34, 3, 144, 4, 48, 136 

(22, 64, 32, 133, 8, 19, 132, 145 

(129,33,3,16,136,68,6,104 

(56, 68, 65, 18, 129, 36, 130, 137 

(48, 64, 3, 144, 40, 129, 14, 133 

(165, 64, 44, 5, 136, 144, 2, 48 

(130, 68, 2, 18, 136, 20, 104, 1 

(13, 74, 40, 133, 98, 193, 18, 148 

(11,65,33, 136, 146,20,152, 100 

(44, 48, 6, 3, 24, 81, 132, 192 

(128,64,32,8, 17,4,2,9 

(3,64,12, 16,33,36,130,136 

(128, 10,68,16,72,36,34, 1 

(42, 64, 6, 16, 9, 37, 132, 129 

(3, 33, 34, 68, 8, 20, 192, 144 

(129, 20, 33, 196, 9, 112, 2, 28 

(3, 65, 164, 16, 192, 34, 12, 72 

(237, 251, 35, 164, 47, 241, 85, 234 

(41, 22, 201, 198, 161, 82, 30, 45 

(6, 64, 32, 16, 136, 3, 140, 129 

(13, 80, 35, 136, 98, 132, 56, 145 

(45,64,6,16,9,160,130,129 

(153,64,32,9, 136,4,2,24 

(25, 65, 160, 67, 136, 4, 130, 48 

(48, 64, 11, 24, 129, 4, 130, 49 

(21, 65, 32, 66, 136, 145, 134, 28 

(20,72,36,194,9,160,6,17 

(25, 33, 194, 20, 136, 164, 34, 68 

(21, 80, 160, 133, 24, 131, 6, 104 

(17, 10,32,84,9,6, 192, 144 

(14, 232, 19, 24, 162, 141, 48, 71 

(33,64,36,17,24,6,130, 136 



(13, 64, 32, 21, 152, 144, 2, 136 

(33,64,160,16,10,4,8,129 

(168, 64, 130, 26, 144, 4, 48, 1 

(14, 72, 36, 144, 138, 48, 192, 1 

(128, 12,32,24, 17,69,2,68 

(130,16,32,196,8,5,18,17 

(18,80,32,135,8,69,130,20 

(13, 65, 32, 20, 136, 144, 2, 80 

(128, 12,3,20,9, 112,34,69 

(28,64,160,136,50,3,56,5 

(132,96,67,16,10,160,40,5 

(164, 33, 136, 16, 73, 5, 2, 224 

(35, 64, 132, 20, 9, 161, 18, 140 

(133, 42, 137, 16, 35, 6, 104, 192 

(137, 72, 98, 129, 19, 4, 34, 152 

(129, 64, 32, 22, 8, 130, 20, 17 

(6,64,32, 129,8, 132,22, 17 

(132, 5, 33, 19, 72, 112, 6, 152 

(3, 64, 130, 144, 136, 4, 48, 41 

(19, 66, 32, 136, 137, 132, 84, 25 

(162, 50, 137, 24, 194, 36, 84, 65 

(129,5,160,84,9,98,18,28 

(5,66,32,144,8,192,6,21 

(44, 96, 132, 19, 136, 81, 130, 21 

(7, 88, 40, 74, 161, 193, 22, 148 

(12, 66, 192, 16, 40, 133, 34, 33 

(128, 166, 172, 72, 107, 50, 85, 113 



(128,64,3, 1,8,4,32,48 

(3,72,32, 144,80,4, 138, 17 

(128,33,40,80,72,4,2, 17 

(132,3,36,20,8,208,34,65 

(25, 64, 32, 129, 152, 4, 2, 137 



R3: M3 ^ C 



R3: U3 ^ C 



16 
15 
17 
16 
13 
15 
15 
16 
21 
17 
15 
19 
17 
18 
17 
15 
17 
16 
18 
17 
17 
15 
22 
20 
18 
10 
14 
13 
16 
15 
18 
16 
38 
28 
14 
21 
16 
14 
17 
16 
19 
17 
19 
20 
16 
25 
15 

16 
12 
15 
17 
14 
14 
17 
15 
17 
18 
16 
17 
18 
20 
19 
14 
14 
19 
15 
19 
21 
19 
15 
20 
23 
16 
27 



10 R3: M3 ^ C 

15 

13 

16 

15 
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12 -30712541418824951 57336 

12 -28706142641542548 57197 

12 -28647676317807325 57054 

12 -26253185098326629 57330 

12 -23870314469588248 57317 

12 -21717939174696339 57321 

12 -21459144106706130 57306 

12 -19374761209641485 57318 

12 -17236351616013344 57315 

12 -15085506227251376 56955 

12 -14842292216336441 56763 

12 -13941487138238279 57270 

12 -9019935303803522 57015 



13 


-15162648248230256 


57309 


13 


-5138871171782917 


57340 


13 


-3841648050655253 


57325 


13 


-1235584894668443 


57338 


13 


202153147298516 


57337 


13 


2399162187229386 


57334 



(41,64,130,144,56,4,10,49) 17 

(136,42,194,80,50,4,224,1) 18 

(128,64,9,10,48,4,18,33) 13 

(5,64,33,152,9,134,130,48) 17 

(3,68,36,20,8,146,192,161) 17 

(130,48,140,26,98,100,76,1) 20 

(128,40,65,16,72,4,2,33) 12 

(29,64,40,137,35,132,18,144) 19 

(28,64,32,136,26,5,130,129) 16 

(14,80,160,133,40,67,146,21) 21 

(3,65,160,80,8,4,34,144) 14 

(131,100,134,24,24,98,196,33) 21 

(42,25,134,69,200,52,162,81) 24 

(128,72,32,24,17,4,2,65) 12 

(22,68,32,144,8,134,192,1) 15 

(14,80,32,73,137,20,146,129) 19 

(128,5,33,17,24,68,2,104) 15 

(12,70,32,193,25,133,18,208) 20 

(12,80,133,18,168,97,6,161) 20 



R3: «4 ^ C 



R3: «3 ^ C 
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